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Abstract. We study the ground state which minimizes a Gross-Pitaevskii 
energy with general non-radial trapping potential, under the unit mass con- 
straint, in the Thomas-Fermi limit where a small parameter e tends to 0. This 
ground state plays an important role in the mathematical treatment of recent 
experiments on the phenomenon of Bose-Einstein condensation, and in the 
study of various types of solutions of nonhomogcneous defocusing nonlinear 
Schrodinger equations. Many of these applications require delicate estimates 
for the behavior of the ground state near the boundary of the condensate, as 
e — > 0, in the vicinity of which the ground state has irregular behavior in the 
form of a steep corner layer. In particular, the role of this layer is important 
in order to detect the presence of vortices in the small density region of the 
condensate, understand the superfluid flow around an obstacle, and also has a 
leading order contribution in the energy. In contrast to previous approaches, 
we utilize a perturbation argument to go beyond the classical Thomas-Fermi 
approximation and accurately approximate the layer by the Hastings-McLeod 
solution of the Painleve-II equation. This settles an open problem (cf. [9] 
pg. 13 or Open Problem 8.1]), answered very recently only for the special case 
of the model harmonic potential |104| . In fact, we even improve upon previ- 
ous results that relied heavily on the radial symmetry of the potential trap. 
Moreover, we show that the ground state has the maximal regularity available, 
namely it remains uniformly bounded in the ^-Holder norm, which is the exact 
Holder regularity of the singular limit profile, as e — > 0. Our study is highly 
motivated by an interesting open problem posed recently by Aftalion, Jerrard, 
and Royo-Letelier |10| , and an open question of Gallo and Pelinovsky 103 , 
concerning the removal of the radial symmetry assumption from the potential 
trap. 
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1. Introduction 

1.1. The problem. This paper is concerned with the analysis of the e — > limiting 
behavior of the Gross-Pitaevskii energy 

Ge{u) = L {} |Vw|2 + i |u|4 + h w[y)H% < ,,y - ' 1 ■ ' 1 

minimized in 

H= \u£ W^ 2 (R 2 ;C) : / W(y)\u\ 2 dy < oo, / \u\ 2 dy=l), (1.2) 

I JR 2 JR 2 J 

where e > is a small parameter and, unless specified otherwise, the potential W 
will satisfy: 

W is nonnegative, W £ C 1 , (1.3) 

and 

there exist constants C> 1, p > 2 such that ^\y\ p < W(y) < C\y\ p if |y| > C, 

(1.4) 

(see also Remark 11.21 below) . It is common to refer to the above problem as the 
minimization of G E under the unit mass constraint. 

Let Ao > infjf2 W(y) be uniquely determined from the relation 

(A -^(y)) + dy = l, (1.5) 

where throughout this paper we will denote / + = max{/, 0}. The choice of the 
value one in the above relation is dictated by the constraint (|1.2I) . see also (|1.13[) 
below. We further assume that the region 

2? = {yeR 2 : W(y) < A } (1.6) 

is a simply connected bounded domain, containing the origin, with smooth boundary 
9£>o, such that 

dW 

— > on 8V 0l 1.7 
on 

where n = n(y) denotes the outward unit normal vector to dT>Q. This last as- 
sumption can be viewed as a non-degeneracy condition. We point out that these 
hypotheses admit physically relevant examples, used to model certain experiments 
(see the next subsection). We stress that the simply connectedness assumption is 
assumed for convenience purposes only (see Remark 11.11 below) , and so is the fact 
that the setting is two-dimensional (see Remark 13. 121 below) . 
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It follows from [186j Prop. 1] (see also [TO], [125], [152]) that the functional G e 
has a unique real valued minimizer r/ £ > in % (all complex valued minimizers are 
of the form rj e e za , where a is a constant). The function r] e satisfies 

- A Ve + ^n E (W(y) + rff) = ^A £ 77 £ , t] £ > in R 2 , 7] s -> as |y| -> oo, (1.8) 

where jrA e is the Lagrange multiplier, which is also necessarily unique. The point 
being that (ll.4[) ensures that minimizing sequences of G e in H cannot have their 
mass escaping at infinity (see also |178j ); in fact the imbedding T~L L 2 (M 2 ,C) is 
compact (see [125] . |208j . or more generally [37] Lemma 3.1]). (One can also ignore 
the mass constraint, and instead minimize the functional G E (u) — ^IMI^Or), w ^ tn 
A > minR2 W, in which case the minimizer would satisfy (|1.8[) with A e = A). 

The real issue is the study of the asymptotic behavior of the minimizer r\ € (or more 
generally of the critical points) of G e as the parameter e tends to zero. Following 
Aftalion and Riviere [4], letting 

A = A - W, (1.9) 
the functional G £ can be rewritten as 

GM = L |v "' ! + i<" 2 - A+ > 2 + i^K+i ( A » - i I <A+)2 * 

(1.10) 

if rj € "H is real valued, where A + = max{A, 0} and A = — min{A, 0}. Let G\{rj) 
denote the first integral above. Since ?7 e clearly minimizes G\ mT-L, by constructing 
a suitable competitor based on V A+, it is easy to see that 

GK%)<^|lne|» (i- 11 ) 
for some constant C > 0, provided e is small (see [8], [TOJ, [130] . and Remark 13.141 
herein). (We remark that the logarithmic term appears because (11.6[) . (11.7[) imply 

that V ^V^A+^j is not square-integrable near dT>o). Hence, for small e > 0, we have 
/ {( n 2 ^A + ) 2 +A- v 2 }dy<Ce 2 \\ne\, (1.12) 

JR 2 

which suggests that i] 2 should be close, in some sense, to A + as e — > 0. Indeed, it 
can be shown that 

i] e -> \/]4+ uniformly in M 2 as e -> 0, (1-13) 

see the references in Subsection 11.31 below. Therefore, loosely speaking, the min- 
imizer rj e develops a steep corner layer along dD 0l as e — > (recall (|1.6[) , (|1.7[> ). 
Note also that V A + is the global minimizer of the "limit" functional 

4e 2?? 2e 2 

among real functions such that ||»7|| z, 2 (h 2 ) = 1 an d Wrj 2 € L 1 (M 2 ). In the context of 
Bose-Einstein condensates, see the following subsection, the function V A + is known 
as the Thomas-Fermi approximation. 

The estimates that are available in the literature for the convergence in (|1.13[) . 
see Subsection II .31 below, fail to encapsulate important information which is often 
required in interesting applications (see Section 11.21 below) . As an illustrative ex- 
ample, let us mention that a lower bound for r\ ei sufficient to imply that the second 
variation d 2 G £ (r] e ) is coercive (this is easy to prove but hard to estimate), does not 



G°(V)= \T^V 4 + 7^W(y)r, 2 \dy, (1.14) 
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seem to be known (in Remark I A. II we will establish a spectral bound for d 2 G £ (rj e ), 
and as a matter of fact one may even calculate sharp constants). Our main goal in 
this paper is to provide crucial details, missing from the known results that describe 
the statement (|1.13j) quantitatively, placing special emphasis on how rj e converges 
to the "singular limit" V A + near DVq, as e — > 0, and proving that it converges in a 
self-similar fashion as conjectured in [9]. Although it is a variational problem, our 
approach will be based more on partial differential equation and functional analysis 
tools. Our treatment is concise and systematic, and can be used to treat in a unified 
manner problems with similar features. 

As will be apparent from one glance in the references of this paper and the 
following discussion, the current interest in the minimizer r] s , and in problems that 
have it on their background, is phenomenal. 

1.2. Motivation for the current work. The motivation for the current work is 
threefold: 

1.2.1. Minimization of a Gross-Pitaevskii energy describing a Bose-Einstein con- 
densate in a potential trap. Among the many experiments on Bose-Einstein con- 
densates (which we abbreviate BEC), one consists in rotating the trap holding the 
atoms in order to observe a supcrfluid behavior: the appearance of quantized vor- 
tices (see the books [9] , j!76] , [177] and the references that follow) . This takes place 
for sufficiently large rotational velocities. On the contrary, at low rotation, no vor- 
tex is detected in the bulk of the condensate. In a BEC, all the atoms occupy the 
lowest energy state so that they can be described by the same complex valued wave 
function. The latter is at the same time the macroscopic quantum wave function 
of the condensate and minimizes a Gross-Pitaevskii type energy. A vortex corre- 
sponds to zeroes of the wave function with phase around it. In two dimensions, the 
Gross-Pitaevskii energy considered in [5] , [10] , [125] , [126] has the form: 

E-M = J |i|Vz;| 2 + ^M 4 + ^W(y)\v\ 2 - fty^ ■ (*«, V«)| dy, v e H, 

(1.15) 

where f2 is the angular velocity, y = (yi,y2) 5 y 1 " = (— y2>yi)) £ > is a small pa- 
rameter that corresponds to the Thomas-Fermi approximation 93, 195 , the trap- 
ping potential W belongs in the class described in the previous subsection, and 
{iv, Vu) = ivVv* — w*Vu. It is clear that rj e is the unique real valued minimizer of 
E e in H. 

For mathematical studies in the case where the condensate has an annular shape, 
we refer to [5] and [53] > whilst for studies in a three-dimensional setting to [T3] , [IS] , 
and |209] . For numerics we refer to the review article [3"U] . 

The density of the condensate is significant in T>q (keep in mind (|1.13p ). which is 
typically a disc or an annulus, and gets exponentially small outside of this domain. 
The case of harmonic trapping potential 

^(yi,y 2 ) = y 2 i+A 2 y2, (1.16) 

for a fixed parameter < A < 1, has been considered in experiments in [15811159] . In 
recent experiments, in which a laser beam is superimposed upon the magnetic trap 
holding the atoms, the trapping potential W is of a different type [1811 11921 1204] : 

W{r) =r 2 + ae- br \ r 2 = y 2 + y 2 , a,b > 0. (1.17) 
(By choosing a, b accordingly, the domain D is either a disc or an annulus). 
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The energy E £ bears a formal resemblance to the well-studied Ginzburg-Landau 
functional 

•h(u) = f |i|Vu| 2 + (\u\ 2 - l) 2 - Qy^ ■ (iu, V«)} dy, 

used to model superconductivity [35] (with il = 0), and superfluidity [13], [188] . 
In their influential monograph [39], Bethuel, Brezis and Helein have developed the 
main tools for studying vortices in "Ginzburg-Landau type" problems. As we have 
already pointed out, the singular behavior of %/ A + near dT>$ induces a cost of order 
| In £■ | in the energy. This causes a mathematical difficulty in detecting vortices 
by energy methods, since any vortex has precisely the same cost (see [39]). In 
other words, the possible presence of vortices will be hidden by the energetic cost 
of the corner layer. This difficulty is common in problems of Ginzburg-Landau 
type when the zero Dirichlet boundary condition is imposed (see for instance [188] ). 
Fortunately, this difficulty can be surpassed in an elegant way by an idea that goes 
back to the work of Lassoued and Mironescu j!48j , and Andre and Shafrir [23] . By 
a remarkable identity, for any v, the energy E £ (v), for any splits into two parts, 
the energy G e {r) e ) of the density profile and a reduced energy of the complex phase 
w = v/r) e : 

E e (v)=G e ( Ve ) + F e (w), (1.18) 

where 

F e (w) = J |VH 2 + ^ (M 2 - I)' - V^Y 1 - ■ iiw, V«>) j dy, (1.19) 

(see also |125j ). In particular, the potential W only appears in G e . In (|1.18[) . 
the term G e {rj £ ) carries the energy of the singular layer near dD$, and thus one 
may detect vortices from the reduced energy F £ by applying the Ginzburg-Landau 
techniques to the energy 




The difficulty will arise in the small density region of the condensate, namely R 2 \2?o, 
where % is small. This kind of splitting of the energy is by now standard in the 
rigorous analysis of functionals such as E e (see also [154] ) . It clearly brings out the 
need for the study of the minimizer rj e of G £ , which is the subject of the current work. 
In particular, as will also be apparent from the discussion in Subsection 11.31 below, 
estimating r\ e near <92?o is essential for adapting the powerful Gamma-convergence 
techniques, developed for J E (see [183] and the references therein), to the study of 
E e (concerning issues of vortices, vortex lines [T3], etc). Obtaining these delicate 
estimates, without imposing any symmetry assumptions on the trapping potential, 
is the main contribution of the present paper. 

1.2.2. Semi-classical states of the defocusing nonlinear Schrddinger equation. Ellip- 
tic problems of the form (|1.8|) arise directly when seeking standing wave solutions 

U (y,t) = e - lAt /^(y), (1.20) 

for the famous nonlinear Schrodinger equation (NLS): 

ie— +e 2 ku-W(y)u±\u\ q - l u = 1 y e R N , t > 0, u : M. N -> C, (1.21) 
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where N > 1 and q > 1 (in the plus sign case, the potential W may vary from that 
described previously). (See also (|6.4[) below). For small e > 0, these standing- wave 
solutions are referred to as semi-classical states. The plus sign in (|1.21|) gives rise to 
the focusing NLS (attractive nonlinearity) , while the minus sign to the defocusing 
NLS (repulsive nonlinearity) which is also known as the Gross-Pitaevskii equation 
(GP). It is also quite common to use the name of the Gross-Pitaevskii equation 
if (|1.21[) has a nonzero potential W, and the name of the nonlinear Schrodinger 
equation if W = 0. Keep in mind that potentials of quadratic growth, as |y| — > oo, 
are the highest order potentials for local well-posedness of (ll.21[) , see |169j . 

At low enough temperature, neglecting the thermal and quantum fluctuations, a 
Bose condensate can be represented by a complex wave function u(y,t) that obeys 
the dynamics of the NLS equation, see the excellent review article [56] and the 
references that follow. In particular, solutions of (|1.8p provide, via (|1.20p . standing 
wave solutions for the GP equation with N = 2, q = 3. Let us mention that the 
minimizer of G £ , considered in the entire space or in a bounded domain with zero 
boundary conditions (as in Remark l3. 131 below), also plays an important role in the 
study of multi-component BECs (see [133] . [155] and the references therein); in the 
dynamics of vortices confined in T> under the flow of the Gross-Pitaevskii equation 
(see [138 ); in the construction of traveling wave solutions with a stationary or 
traveling vortex ring to the Gross-Pitaevskii equation (see [180 , [202], |203] ); in the 
study of excited states of the GP equation (see |119j . [174] ). and in Bose-Einstein 
condensates with weak localized impurities (see |100 ). 

In considering typical BEC experiments and in exploring the unprecedented con- 
trol of the condensates through magnetic and optical "knobs" , a mean-field theory 
is applied to reduce the quantum many-atom description to a scalar nonlinear 
Schrodinger equation (see |152j ). The NLS equation is a well established model in 
optical and plasma physics as well as in fluid mechanics, where it emerges out of 
entirely different physical considerations [Tl 1193] . In particular, for instance in op- 
tics, it emerges due to the so-called Kerr effect, where the material refractive index 
depends linearly on the intensity of incident light. The widespread use of the NLS 
equation stems from the fact that it describes, to the lowest order, the nonlinear 
dynamics of envelope waves. 

Ground state solutions of the NLS are standing wave solutions, of the form (II -20[) . 
such that rj is positive, r\ 6 W 1,2 (R N ), and satisfies 

e 2 A?y - (W(y) - A) r) ± M 9-1 ?? = in R N , rj ->■ as |y| -> oo. (1.22) 

The condition u € W l ' 2 (M. N ) is required to obtain solutions with physical interest. 
(Sometimes we will refer to positive solutions rj of (|1.22|) as ground states of (|1.2ip ). 
In the subcritical case where 1 < q < if TV > 3, q > 1 if TV = 1,2, ground states 
of the defocusing equation (|1.22p _ correspond to global minimizers of G- (these are 
nontrivial if e > is sufficiently small |125j . see also [HJ Example 5.11] and [821 
Lemma 2.1]), while ground states of the focusing equation (|1.22[) + correspond to 
mountain passes of G+, where 

G±(V)=[ \ £ 4\^\ 2 + (W(y)-X)^T l -^}dy, (1.23) 
Jr n 12 2 q + 1 J 

in W 1 - 2 (R N ) with J RN W(y)r] 2 dy < oo (see for instance [H], [T75] V 

In the focusing case, following the pioneering work of Floer and Weinstein [98], 
there have been enormous investigations on spike layer solutions for problem p.22p + . 
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for small £ > 0, typically under the conditions A < inf R w W(y) < liminf | y |->-oo W(y), 
and 1 < q < $±§ if N > 3, q > 1 if N = 1, 2. The "critical" case, where 
inf R jv W(y) = A, has also received attention, see [17] ■ Actually, the latter case is 
related to the discussion following Definition [1] below. We refer the interested reader 
to [2D], [75], and the references therein. 

In the defocusing case (ll.22[) _, assuming that W satisfies the assumptions of 
the previous subsection with A in place of Ao, we see that the corresponding limit 
algebraic equation, obtained by formally letting e — in (11.22p _ , has the compactly 
supported continuous solution 

Vo = [(A - W)+] . (1.24) 

Obviously, there is also the solution — 770 , and the trivial one. In fact, if q > 2 then 
the "singular limit" 770 is merely Holder continuous with exponent l/(q — 1) (recall 
(ll.6[) . (|1.7[) ). In particular, it is easy to see that if q > 3 then V770 is not square- 
integrable near <92?o- If g = 2 and 77 > 0, then (|1.22l) _ becomes the well known 
scalar logistic equation [55], and 770 is Lipschitz continuous. If 1 < q < 2, then 770 is 
at least differcntiable. In the language of bifurcation theory [128] . the solution set 
of the corresponding limit algebraic equation to (|1.22p _ undergoes a supercritical 
pitchfork bifurcation at dT>o if q > 2; a transcritical bifurcation if q — 2; a subcritical 
pitchfork bifurcation if 1 < q < 2. Notice also that 770 is asymptotically stable (as an 
approximate equilibrium) with respect to the parabolic dynamics that correspond 
to (|1.22[) . The question whether 770 perturbs, for small e, to a solution n e of (|1.22[) _ . 
and keeping track of the maximal regularity available (to be defined in a moment), 
is a source of current mathematical interest. Note that, if q > 2, such a solution 
would have a corner layer along &Dq. The bifurcation that occurs at 8T>q takes 
the problem off from the classical setting, where the roots of the corresponding 
algebraic equation (for fixed y) are non- intersecting (see |77|). 
The following definition is adapted from [52/. 

Definition 1. Let a € (0,1] be the largest number such that ||77o||c ' Q (R N ) < °°; 
we say that a solution n E has maximal Holder regularity if ||??e||c ' O! (M JV ) remains 
uniformly bounded as e — > 0. 

If N = 1, problem (|1.22[) _ can be rewritten as a homoclinic connection problem 
for a 3-dimensional slow-fast system of ordinary differential equations [131] , At 
the points that correspond to dT>o we have loss of normal hyperbolicity of the 
corresponding slow manifold, due to a pitchfork or transcritical bifurcation. This 
fact prohibits the use of standard geometric singular perturbation theory |92| 1131] 
in order to deduce the persistence of the "singular homoclinic orbit" 7/0, for small 
£ > 0. At present, much work in geometric singular perturbation theory deals with 
such situations, often using the "blowing-up" construction [145] , This approach 
has been successfully applied recently in 190^ in a heteroclinic connection problem 
for a 4-dimensional slow-fast Hamiltonian system, sharing similar features with 
(|1.22|) _ with q = 3, arising from the study of crystalline grain boundaries (see also 
[18 ] l97l T191] ). Let us mention that there is an abundance of non- hyperbolic points 
in applications, see for instance [145] . that have been successfully treated using this 
approach. In particular, singularly perturbed one-dimensional second-order elliptic 
systems, involving loss of normal hyperbolicity, arise in the study of the Dafermos 
regularization for singular shocks [66], [189] . Although elegant, the arguments of 
this approach are intrinsically one-dimensional. 
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Elliptic systems where the singular limit has merely Holder or Lipschitz regular- 
ity typically describe phase separation or spatial segregation, and have attracted a 
lot of current mathematical research, see [52j |62j 199 and the references therein. 
These type of problems have been tackled in the latter references either by weak 
convergence arguments, yielding weak convergence in the Sobolev space W 1,2 and 
strong in L 2 (see also 69 ), or the method of upper and lower solutions, yield- 
ing uniform convergence (see also [IB], [123] . [127] ). as e — > 0. The question of 
maximal regularity of the convergence to the singular limit profile is then addressed 
using a delicate analysis, based on monotonicity properties, blow-up techniques and 
Liouville-type theorems (see also Remark 14.11 below) . To the best of our knowledge, 
for these systems, the maximal regularity property has only been proven in one- 
dimensional cases, see [36] . In the case at hand, since Vryo is not square-integrable 
near dT>o if q > 3, it is not clear how to use standard weak convergence argu- 
ments. Furthermore, it seems to be hard to construct a good pair of upper and 
lower solutions (especially) near 8T>q. A motivation for the current work is to show 
that the perturbation approach to such problems, we initiated in [134U135] . can be 

adapted to treat problem (II . 22[) with general potential. We emphasize that the 

perturbation method seems to be the only one available at the moment that yields the 
maximal regularity property in higher- dimensional singular perturbation problems. 
(See the main theorem of [135] . and Corollary 14 . 31 herein) . Among other advantages 
of the perturbation approach is that it provides finer estimates which, as we already 
stated, imply the maximal regularity property, it be applied in the study of non- 
minimizing solutions in systems of equations which in general lack the maximum 
principle (see assumption (B4) in [45]), and in supercritical problems that cannot 
be treated variationally (see for example Remark IB.5|) . Most importantly, for the 
problem at hand, it applies equally well without the assumption of radial symmetry 
on the equation, see the discussion in the following subsection. 

1.2.3. Applications to related problems. A strong motivation behind the current 
work is the possibility of adapting our approach to treat other interesting problems 
with similar features. 

Variational problems of a similar type to (|1.15p . with spatially varying coeffi- 
cients, have also been introduced to model vortex pinning due to material inho- 
mogeneities in a superconductor (see the references in [12], [83] Sec. 6] and [1831 
Subsec. 14.1.4, 14.3.5]). In [22], the authors considered the case where the corre- 
sponding pinning potential W(y) is nonnegative but vanishes at a finite number of 
points (see also [194j ). Minimizers such that |u e | — > vVV in W 1 ' 2 , as e — > 0, were 
analyzed by variational methods. In their result, it was important that vVV is in 
W 12 , which is not the case here (also recall the previous discussion concerning weak 
convergence methods). In that context, the minimizer has corner layer behavior at 
points rather than curves. It would be of interest to study this type of problems 
from the perturbation viewpoint of the present paper which, in particular, does not 
require that the singular limit belongs in W 1 ' 2 . 

A singularly perturbed elliptic equation of the form (|1.50|) below, considered in 
a bounded domain with zero boundary conditions, where the corresponding limit 
algebraic equation admits a fold bifurcation at the boundary of the domain, appears 
in the proof of the Lazer-Mckenna conjecture for a superlinear elliptic problem of 
Ambrosetti-Prodi type, see [71] . [73j . In that case, the corresponding minimizer 
(without the mass constraint) has a steep corner layer along the boundary of the 
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domain, see also the old paper |197j . This situation is qualitatively similar to the 
problem considered here. In [71] . by adapting variational techniques from [70], the 
behavior of the minimizer, as e — > 0, was estimated in compact sets away from the 
boundary of the domain. We believe that, employing the perturbation techniques of 
the present paper, one can obtain fine estimates for the minimizer all the way up to 
the boundary. In turn, these could potentially lead to the construction of new type 
of solutions on top of the minimizer, for example solutions having small peaks near 
the boundary, or bifurcating from symmetry (as in |134[ 1135] ). Another reason for 
developing perturbation arguments for these problems is that unstable solutions are 
hard to find or describe accurately through purely variational methods, see Remark 
IB. 51 for more details. 

I. 3. Known results. In this subsection we gather some known properties of the 
real valued minimizer rj e of G e in H, under the assumptions described in Subsection 

II. 11 These have been studied in various contexts (see for instance [8j [TO] EU I125( 
I163j ). As we will see in this paper, some are actually far from optimal. 

The corresponding Lagrange multiplier satisfies, for small e > 0, the estimate: 

|A S — AqI < C7j lne|*e. (1.25) 
The real valued minimizer r\ e satisfies the following estimates: 

< r] e < VA+ + Csi inR 2 , (1.26) 

Ve(y) < C*£*exp{-c£-idist(y,<9P )} in M 2 \2? , (1.27) 
where y dist(y, dT>o) denotes the Euclidean distance in M 2 from y to dT>o, 

|ry £ - VA + \ < Ce^VA+ if y <= V and dist(y, <9X> ) > e% , (1.28) 

and 

HVTfelUo.^ < Ce~\ (1.29) 
for some constants c, C > 0, if e is sufficiently small. Relation (jl.25l) follows from 
([TTTT]) . see [10], [125]. Relations ([L~26]t . (fl~27f have been shown in [125] (see also 
[163( Lemma B.l]) by constructing a suitable upper-solution of (|1.8j) . and using 
the uniqueness of positive solutions of the latter equation in bounded domains with 
zero boundary conditions [23] . Relation (|1.28[) can be traced back to [5] , and follows 
from lower- and upper-solution arguments in equation (jl.8l) based on [23 (also keep 
in mind Remark 3 in [8]). Note that, in particular, estimates (|1.26[) - (|1.28j) yield 
(|1.13p . Lastly, estimate (|1.29j) on the gradient follows from the equation and a 
Gagliardo-Nirenberg type inequality as in [38] (see Lemma 14.11 herein) . 

In the special case where the potential trap W is additionally assumed to be 
radially symmetric, it follows from its uniqueness that the real valued minimizer 
T] e > of G e in T-L is also radially symmetric. In particular, if T>q is a ball of radius 
R, it has been shown recently in [10] that 

i£(r)<0 m(R-S ,R + S ), (1.30) 

for some small constant Sq > 0, if e is small. The radial symmetry was used heavily 
by the authors of [lOj in order to establish (|1.30[) , using a maximum principle due to 
Berestycki, Nirenberg, and Varadhan [3S] (see also the discussion following (|3.29l) 
herein) together with an intersection-comparison type of argument, mostly taking 
advantage of (|1.7j) . The importance of the positivity of % and the radial symmetry 
of W in deriving (|1.30j) can be naively seen from the following consideration. If 
W'(r) > for all r > 0, using that n e > and the method of moving planes [1081, 
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we can infer that rf e (r) < for all r > (see also Proposition 2.1 in |152) for an 
approach via a radially-symmetric rearrangement argument which takes advantage 
of the minimizing character of rj e ). We point out that it is not clear how the 
aforementioned arguments apply in the case where W is radially symmetric but the 
set {W — Ao < 0} is an annulus, considered in jS], [13] • Relation (|1.30l) was used in 
an essential way in |10j for estimating uniformly the auxiliary function 

f £ (r)=Ur)/Ve(r), where &(r) = / sii(s)d Sl (1.31) 

J r 

near the circle |y| = R. The function f £ plays a crucial role in the study of the 
functional E e , see [TU] and Section [3] below. (Actually, estimate (|1.30l) was needed 
in a region of the form (R — S, R + Ce^) for some constants S, C > 0). Making use 
of the previously mentioned estimates on f £ near the boundary of T>o, and of some 
new estimates away from T>$ > the authors of [10] proved that if the angular velocity 
is below a critical speed 

fl c = uio\ hie| — u)% In | lne|, 

for some constants loq, uj\ > 0, and e is sufficiently small, then the rotation has 
absolutely no effect on the minimizer. In other words, all minimizers v £ of E £ in 
T-L are of the form v £ — r\ £ e %a in R 2 , where a is a constant. In particular, at low 
velocity, there are no vortices in the condensate. We remark that this last assertion 
was previously known to hold true only in the bulk of the condensate, see |125j . It 
was left as an interesting open problem in |10j to see to what extent their analysis 
continues to hold if the assumption of radial symmetry on W is dropped (see also 
Open Problem 8.1 in [S], and the open questions in the presentation [103] ). Our 
results on the minimizer r] £ , which hold without any symmetry assumption on W, 
may represent a major step in the answering of this question. 

Let us close this subsection by mentioning that the case where the potential is 
homogeneous of some order s > 0, i.e., W(Xy) = X s W(y) for all A > 0, y 6 R 2 (see 
(|1.16p for an example with s = 2), and locally Holder continuous has been studied 
in the work of E. H. Lieb and his collaborators in |152j . [153] . By employing scaling 
and variational arguments, it was shown in the latter references that, as e — > 0, the 
minimizer % converges to y/ A + in the strong i 2 (R 2 ) sense. In the special case of 
the model harmonic potential, described by (|1.16[) with A = 1, a complete analysis 
has been carried in [104] (see the discussion following the statement of our main 
theorem for more details). By generalizing the divergence- free WKB method, a 
uniformly valid approximation for the condensate density of an ultra-cold Bose gas 
confined in a harmonic trap, that extends into the classically forbidden region (near 
and outside of &Dq), has been established very recently in |182j . 

1.4. Statement of the main result. In order to state our main result, we need 
some definitions: 
Let 

Oe(y) = A e - W(y), yeR 2 . (1.32) 

By virtue of (ll.6[) , (ll.7[) , (ll.25|) , and the implicit function theorem [115] , the domain 
T>q perturbs smoothly, for e > small, to a simply connected domain T> £ such that 

a £ >0 in V £ , a £ < in R 2 \2? e , < -c on &D £ , (1.33) 
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for some constant c > 0, where v e = v e (y) denotes the outward unit normal to dT> £ . 
Let T £ be the simple, smooth closed curve defined by dV e , and let 4 — \T £ \ denote 
its total length. Note that r e is inside of an 0(\ lne| 2 e)-tubular neighborhood 
of dT> . We consider the natural parametrization 7 e = r ) £ {9) of T E with positive 
orientation, where 9 denotes an arc length parameter measured from a chosen point 
of T e . Slightly abusing notation, we let v s {ff) denote the outward unit normal to 
r e (as in (|1.33p ). Points y that are 5o-close to T e , for sufficiently small 5q > 
(independent of small e), can be represented in the form 

y = -y e (e) + tv e {0), \t\<6a, 0e[O,4), (1.34) 

where the map y H> (t, 9) is a local diffeomorphism (see )109| Sec. 14.6]). Note that 
t < in T> s . Abusing notation, frequently we will denote points y near T e plainly 
by their image (t,9) under the above mapping. From (|1.32|) . (|1 .33|) . we have 

-(a e ) t (0,9)=Wt(0,e)>c, 6 [0,4), (1.35) 

for some constant c > and small e. We define 

f3 e (0) = (-a t (0,9))i > 0, 9e [0,4)- (1-36) 

It might be useful to point out that for the harmonic potential, described in (|1.16j) . 
one can derive explicitly that 



[/3 £ (0)] 3 -2 v ^7A^A-2 C os2 +sin 2 (j-O) , G [0,4). 

We will also make use of the Hastings- McLeod solution |113j of the Painleve-II 
equation [99 , namely the unique solution V of the boundary value problem: 

v xx —v{v 2 +x) =0, x & M, v(x) — \J — x — >• 0, x — » — 00; v(x) — > 0, x — > 00. (1.37) 

It is useful, in relation with (ll.30|) . to point out here that V x < 0, x € R. The 
importance of the Hastings-Mcleod solution is that it will "lead" the minimizer rj e 
across dT> . 

We can now state our main result: 

Theorem 1.1. If e > is sufficiently small, the unique real valued minimizer of 
G £ in H satisfies 

% (y) = sif3 e (0)V (p e (6)-t) + 0(e + \t\i) (1.38) 



£3 



uniformly in {-d <t<0, 9 € [0, 4)}; 

r, e (y) - e*[3 e (0)V (^W4) + ©(ee"^ 1 ') (1.39) 

uniformly in {0 < t < d, 9 € [0,4)}; where c, d > are some small constants. 
Given D > 0, if e > is sufficiently small, we have 

Ve(y)-V^-W(y)=o(e 2 \t\-^ (1.40) 

uniformly in {-d<t<-De%, 6»e[0,4)}, 

% (y)-^/\ e ~W(y)=0(s 2 ) (1.41) 
uniformly in T> e \ {—d < t < 0, 9 S [0,4)}; an d 

< rj e (y) < Ce^ exp{-c£ - § dist(y, dT> )} (1.42) 
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in M. 2 \D e , for some constants c, C > independent of e, where O(-) denotes Lan- 
dau's symbol (see Section^ for the precise definition). 
The potential of the associated linearized operator 

L e (<p) = e 2 Ap - (3n 2 (y) + W(y) - A £ ) ip, (1.43) 

satisfies 

( cei+c\t\, if \t\<5, 
3 V 2 e (y) + W{y) -\ e >l (1.44) 
[ c + c|y| p , otherwise. 

The Lagrange multiplier e~ 2 A £ satisfies 

A e -A -O(|lne|e 2 ), (1.45) 

while the energy of n e satisfies 



gm L {A)dy ) 6 + u[ J l ,,, -"l- c ' 1| ' il - (L4,,i 

as e —> 0. 

The main highlight of our result is that we rigorously prove that, close to dT>o, 
the minimizer r\ e behaves like (|1.38l) . We emphasize that the rigorous derivation of 
the Paileve-II equation from (|1.8I) was an unsettled open problem, see [S] pg. 13] 
(see also the discussion below on the recent paper |104) ). In turn, as was noted 
in Section 8.1.3 of the book [S], this information is required as a stepping stone 
towards the open problem mentioned in Subsection 1 1.31 (see also Open Problem 8.1 
in |10j). in order to obtain a lower bound for r\ e in R 2 \2? - I n the current paper, 
under the additional assumption that Wu is strictly positive on dT>o, we contribute 
towards this direction by obtaining a sharp lower bound in the strip-like domain of 
R 2 \r> described by dist(y, dV a ) < si (see Remark S3] below). 

We believe that our result opens the road for the rigorous description of the 
Painleve region in recent experiments on three-dimensional Bose-Einstein conden- 
sates where a laser beam, modeled by a cylinder along the z direction, is translated 
in the x direction along the condensate (see [61(7]^, and to understand the superfluid 
flow around an obstacle (see [9] and the references therein). 

As we have already mentioned, our proof is based on perturbation arguments 
rather than variational ones or the method of upper and lower solutions (as was 
hoped for in [5J Sec. 8.3.1]). It relies mainly on the existence and asymptotic 
stability of the Hastings-McLeod solution (in the usual sense) in order to construct 
a solution u e of problem (jl.8p . "close" to V A + , for small e > (with the Lagrange 
multiplier A e treated as a known coefficient). Then, using the fact that the latter 
problem has a unique solution which follows from ideas of Brezis and Oswald |43] 
(see Remark 13.121 below) , namely r/ e , we infer that u e = n e . Actually, with the 
obvious modifications, an analogous result holds true for the minimizer of G e in 
arbitrary dimensions. Furthermore, our method of proof can be adapted to treat 
the case where G e {u) — ■|||m|| 2 2( R 2- ) is minimized in Wq' 2 (D), where ED is a bounded 
annular domain such that W = A on the outer or inner part of its boundary, as in 
[5], [13], see Remark 13.131 below. For further generalizations we refer to Remarks 
EDO below. 

It follows from the above estimates that the convergence of r\ e to V A + is the 
most regular possible (see Corollary 14.31 below) . If we assume that W u > c > on 
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r, then we will show in Proposition 14. II that bound (|1 .39|) can be replaced by 
r, e (y) = e^MV (&(0)4) 

Estimate (|1.45p improves (|1.25l) . and was previously established in |125j in the 
special case of the model harmonic potential (|1.16[) by exploiting a scaling property 
of the corresponding equation (ll.8[) . see Remark 13.151 below . which is not available 
under our general assumptions. The above theorem has some other interesting 
consequences, which will be explored in Sections 01 [5] We can prove an analogous 
monotonicity property to (|1.30[) for r\ e near dT>o without the simplifying assumption 
of radial symmetry on the potential W. We will see that estimates (|1.28[) and (|1.29[) 
are actually far from optimal. In addition, restricting ourselves to the radially 
symmetric case with T>q a ball, we can improve and sharpen the new estimates of 
[TU] for the important auxiliary function f e , as described in (|1.31[) . In fact, we believe 
that the estimates of Theorem 11.11 can be utilized in estimating the function f e in 
the non-radial case, through equation (|5.ip below, which may ultimately lead to 
the resolution of the open problem raised in [10] (recall the discussion in Subsection 

El. 

To further illustrate the importance of our result, we emphasize that its method 
of proof can be adapted to produce similar results for semiclassical standing wave 
solutions of the defocusing NLS (|1.2ip _. with N > 1 and q > 2, assuming that W 
has the features described Subsection 11.11 with A in place of Ao and T>o a domain 
with (TV — l)-dimensional boundary (see also Remark IB . 1 1 below ) . (Recall that the 
minimizer r\ e solves (|1.8|) . and the discussion in the second part of Subsection II. 2|) . 
We also emphasize that, when dealing directly with (ll.22[) _ . our approach does not 
make any use of the techniques of Brezis and Oswald [43 , which were needed in 
previous approaches for establishing (jl.26l) . In fact, our approach may produce sign 
changing solutions of (ll.22p ^, satisfying estimates analogous to those of Theorem 
11.11 (sec Remark ll. II and Section [6] below). 

A rigorous connection between semiclassical ground states of the defocusing non- 
linear Schrodinger equation (ll.21[) _. in one space dimension, and solutions of the 
Painleve-II equation (|1.37[) has been established recently in |190l I191j , for a related 
Hamiltonian system. We refer to the physical works [26], [143] . |156] . [176j for for- 
mal expansions in one space dimension or radially symmetric cases, and to [205] 
for higher dimensions (see also [162 ). In the case of the model harmonic potential 
W(y) — |y| 2 , y € R , N > 1, A = 1, and q — 3, the problem of uniform asymp- 
totic approximations of the ground state of the defocusing NLS (|1.21[) _ . using the 
Hastings-McLeod solution of the Painleve-II equation, has been established on a 
rigorous level very recently in [104] . However, the approach of |104j relies crucially 
on the specific form of the model harmonic potential, which allows for a suitable 
global change of independent variables in equation (|1.22p „ (see Remark 13.171 be- 
low). The real delicacy of our result is not that it successfully connects r\ e with V, 
but that we do so in a way that yields fine estimates, as can already be seen from 
(|1.46p , and the proof of (|1.45p (see also the detailed estimates of |104j in the case 
of the model harmonic potential). The optimality of estimates (jl.38p . (|1.40p is also 
suggested by Remark 13.101 below. 

The analog of relation (|1.44l) may be used to study the spectrum of the lineariza- 
tion of the Gross-Pitaevskii equation (|1.2ip _ at the corresponding standing wave 
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solution (|1.20p . which is defined by the eigenvalue problem, in L 2 (R N ) x L 2 (R N ), 
-e 2 A<j>+(W -\ + q\r} e \ q - l )<j) = -/j,ip 

(1.47) 

-e 2 A4>+{W - X+Ij^- 1 )^ = lief). 

The above eigenvalue problem determines the spectral stability of the standing 
wave, with respect to the time evolution of the GP equation, and gives preliminary 
information for nonlinear analysis of orbital stability or more generally about the 
flow of (TPTj) ^ in a neighborhood of the standing wave (see [2S], Q3U, [33D], QH])- 
More complex phenomena, such as those of pinned vortices (dark solitons) on top 
of the ground state [56], |173j . and the construction of traveling wave solutions 
with a stationary or traveling vortex ring to the GP equation [202] . [203] . can 
also be understood from the analysis of (|1.47l) . In particular, relation (|1.44l) plays 
an important role in the construction of excited states for (|1.22P - (see [174] and 
Section [S] herein) . It seems that it was previously known only in the special case 
of the model harmonic potential (see |104j ). In fact, even under the assumption of 
general radial symmetry, it does not follow from the recent estimates of [10] . 
Observe that the equation in (jl.22|) -. with q = 3, is equivalent, for e > 0, to 

Au+ (X-W(sy))u-u 3 , yeR N , (1.48) 

which when setting e — 0, and re-scaling appropriately (we assume that A > W(0)), 
becomes the well known Allen-Cahn equation 

Av - v{v 2 - 1) = 0, ye R N . (1.49) 

The above problem has received an enormous attention, see for instance [81] and 
the references therein. It seems plausible that our result can be combined with 
existing ones for the Allen-Cahn equation, and produce new interesting solutions 
for the Gross-Pitaevskii equation (|1.22[) _. We will elaborate more on this in Section 
[6] below. 

Remark 1.1. Analogous assertions to those of Theorem \1.1\ hold true if T>q is 
assumed to be the union of finitely many bounded smooth domains. In the latter 
case, one may construct sign changing solutions of (|1. 2211 - . whose absolute value 
converges uniformly to V A + , as e — >• 0. The assumption thatT>o is simply connected 
plays an important role only in Section [5| below. The degenerate case where some 
connected components ofT>o "touch" is left as an open problem is Section^ 

Remark 1.2. With only minor modifications in the proof all the assertions of 
Theorem \l.l\ remain true if \l-4ty is replaced by ^ < W(y) — Xq < C\y\ p , \y\ > C, 
for somep > (the second branch in jl-44^ would have to be replaced by a constant). 
The fact that p > 2 in \l-4\) plays an important role only in Proposition I5.il below. 

Remark 1.3. Let us keep in mind that \t\ = disi(y,92?o) + 0{\ In e| ^e), wherever 
defined, and a-posteriori \t\ = dist(y,dT>o) + 0(\ lne|e 2 ), as e — > (recall il.25\) . 

Remark 1.4. An analogous result continuous to hold for the singularly perturbed 
elliptic equation 

e 2 Au = F{u,y), (1.50) 
considered in the entire space M. N or in a bounded domain, where F is such that the 
zero set of F(u, y) = undergoes a supercritical pitchfork bifurcation as the variable 
y crosses some surface (such conditions on F can be found in [190] j. In order to 
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bring out clearly the underline ideas we refrain from any such generalization. In 
il.39\) the convergence is exponential because, in the case at hand, zero is a solution 
of F = 0. However, there are no delay phenomena present, as y crosses dT>o, in 
contrast to the first order equation su — F(u, t), as t crosses the point corresponding 
to dV (see [55] ). 

1.5. Outline of the proof and structure of the paper. The proof of Theo- 
rem 11.11 consists in showing that there exists a genuine solution of (|1.8j) "near" a 
suitably constructed smooth approximate solution, which in turn is "near" V A + , 
provided the parameter e is small enough. We emphasize that we will consider 
the Lagrange multiplier e~ 2 X £ , corresponding to the minimization of G e in TL, as a 
known coefficient in fl.8\) which, as we have already remarked, is known to satisfy 
U.25\) . Then, by uniqueness (see also Remark 13.121 below), we will conclude that 
the obtained solution is actually the minimizer n £ . Actually, we will prefer to work 
with the equivalent (for e > 0) problem in stretched variables y = e _3 y, see (|3.ip 
below. The main steps of the proof are the following: 

(i) Firstly, we construct a "good" smooth approximate solution for the (stretched) 
problem which we call u ap - The function u ap is carefully built, along the 
lines set in [135] . throughout Subsections I3.HT3731 in the following steps: 
Starting from the Hastings-McLeod solution V, described above, we con- 
struct an inner approximation Ui n that is valid only in a tubular neighbor- 
hood of the (stretched) curve e~3T £ . Then, we obtain the desired (global) 
approximate solution u ap by patching u; m with a subtle modification of the 



outer approximation y a £ (e^y) + . We emphasize that the use of this latter 
modification seems to be a key point in the whole construction. In fact, 
the matching of the inner approximation with the outer one is the major 
difficulty in the current problem, mainly due to the algebraic decay of V to 
yj—x as x — > — oo (see (|3.18l) below). For more details on this point, see 
the discussion in Subsection 13.21 and Remark 13.81 below. This difficulty was 
not present in |104j . where the case of the model harmonic potential was 
considered, since in that case the special form of the equation allowed for 
the inner solution to be used globally, and thus no matching was necessary, 
(ii) Next, in Subsection 13.41 we study the linearized operator 



about the approximate solution u ap , and invert it in carefully chosen weighted 
spaces. We exploit a recent observation, due to Gallo and Pelinovsky |1041 
Lemma 2.2], which says that the potential 3V 2 + x of the linearization of 
(|1.37|) about the Hastings-McLeod solution is bounded below by some pos- 
itive constant. In turn, we will show that this latter property implies that 
the same assertion holds true for the potential of — C £ , if e > is small. 
In particular, for small e, the operator C £ is invertible. We point out that 
the previously mentioned lemma of |104) is of technical nature and can be 
bypassed at the expense of a more involved, but rather standard, analysis 
(see Remark 13.71 and Appendix IA1 below). Our choice of weighted spaces, 
a variant of those considered in jl71) . is mainly motivated from the error 
term ([535} below. 
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(iii) Finally, in Subsection 13.51 we look for a genuine solution of the stretched 
problem (13. lj) in the form 

u e = u ap + <p, 

where tp is a correction. At this stage, we show that we can rephrase the 
problem as a fixed point problem for ip, which can easily be solved, if e > 
is sufficiently small, using the fixed point theorem for contraction mappings. 
Then, in Subsection 13.61 taking advantage of the recent uniqueness result 
of |125j for problem (|3.ip (see also Remark 13.121 below) . we infer that the 
unique real valued minimizer of G e in H satisfies ij £ (y) = « e (e"3y), y G M 2 . 
From this property, and the estimates derived from the construction of u £ , 
we can deduce the validity of Theorem 11.11 In particular, the proof of 
estimate (|1.45l) builds on (|1.25l) . and uses the fact that ||?7£||l 2 (r 2 ) = 1- 
The outline of the paper is the following: In Section [5] we will introduce notation 
and standard concepts that we will use throughout the paper. In Section [3] we will 
present the proof of our main result, as outlined above. In Section^ as a byproduct 
of our construction, we will establish an analogous monotonicity property to (| 1 . 30[) 
for the general case (with improvements), show that rj E has the maximal Holder 
regularity available, improve bound (jl.29j) . and generalize and considerably improve 
estimate (|1.28j) . In Section [5j assuming that the potential W is radially symmetric 
with 2?o a ball, we will mainly rely on the results of Section |4] to refine and improve 
the recent estimates of [10] for the important auxiliary function f e in (jl.3ip . In 
Section [6] we will present some interesting open problems that are related to the 
current study. We will close the paper with two appendixes. In Appendix [X] we 
will reprove our main result concerning the linearized operator C s based on the 
asymptotic stability of the Hastings-McLeod solution (in the usual sense), rather 
than making use of the recent lemma of [104] which may not hold in other problems. 
Lastly, in Appendix [B] we will provide a new more flexible and simple proof of the 
existence, and related properties, of the Hastings-McLeod solution through the 
study of problem (|3.107j) below. The results of the latter study seem to be new and 
to have interesting applications (see Remark 13.131 below). 

2. Notation 

In the sequel, we will often suppress the obvious dependence on s of various 
functions and quantities. Furthermore by c/C we will denote small/large generic 
constants, independent of e, whose value will change from line to line. The value of 
e will constantly decrease so that all previous relations hold. The Landau symbol 
0(1), e — > 0, will denote quantities that remain uniformly bounded as e — > 0, 
whereas o(l) will denote quantities that approach zero as e — > 0. By C < a or 
equivalently C k+a , W k > p , k,p G Z U {oo}, a G [0,1), we will denote the usual 
Holder, and Sobolev spaces respectively (see for instance [109] ). Frequently, we will 
denote the minus case of (| 1 . 22|) below plainly by (jl.22|) „. e.t.c. 

We will identify the tubular neighborhood Bs(T E ) = {y G R 2 : dist(y, r e ) < 8}, 
where 8 < 8 (recall ([L34]) ). with Q 5 (T e ) = {\t\ < 8, 6 G [0,4)}, denoted simply by 
{\t\ < 8}, and y in Bs(T e ) by the corresponding pair (t, 9) G Sl,5 (r e ) as determined 
via (fT3% 



3. Proof of the main result 
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3.1. Setup near the Curve. In this subsection, suitably blowing up (|1.8p around 
the curve T e , we will construct an inner approximation which is valid only near the 
curve, for small e > 0, whose profile normal to the curve will be that of a (scaled) 
Hastings-McLeod solution (see also [205 ). To this end, we will follow the general 
lines set in [75] which dealt with the focusing case (|1.22p + . 

Formally neglecting the term e 2 Au in (ll.8|) . we get the outer approximation 
y/a e (y) + . However, the Laplacian of the latter is not even bounded on the boundary 
of T> e . Hence, the outer approximation fails in the vicinity of dT>Q. An inner 
approximation is thus needed, playing the role of a "bridge" as y crosses that 
boundary. 

In the coordinates (t, 6), near T e , the metric can be parameterized as 

g tfi = dt 2 + (l + kt) 2 d6 2 , 
and the Laplacian operator becomes 

d 2 1 d 2 k d k't d 



dt 2 (l + kt) 2 d6 2 l + ktdt (l + kt) 3 d9' 

where k e (8) is the curvature of r e (see for instance [TTJ, [§5])- Note that k £ and its 
derivatives depend smoothly on e > 0. 

Remark 3.1. (t,^/ £ (9)) are known in the literature as Fermi coordinates, see for 
instance |140j . and are frequently employed in the study of perturbation problems 
involving solutions concentrating on manifolds (see for example |78| ). Interestingly 
enough, they owe their name to the physicist E. Fermi in the title of the current 
paper! 

In stretched variables 

V = £~*y, 

problem (| 1 .8[) becomes 



Am - s'iu (u 2 - a £ (e%y)) = 0, u > in M 2 , lim u(y) = 0, (3.1) 
\ / \y\->-oo 

where a e was defined in (|1.32[) . 
In the sequel we will denote 

V £ = s'iv e and f £ = e^Y e . (3.2) 

Let 

{s,z)=e-i{t,6) (3.3) 
be natural stretched coordinates associated to the curve r e , now defined for 

s e (-5 £-% , 5 e-%Y ze 0,£-?4)- (3.4) 
In the coordinates (s, z), near F e , the metric can be written as 

9s,z = ds 2 + (1 + e^ksfdz 2 , (3.5) 
and the Laplacian for u expressed in these coordinates becomes 

Au = u zz + u ss + Bi(u), (3.6) 
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where 
B 1 (u) = - 



1 



ezk(s3 z) 



sk'(e 3 z) 



1 + £3 k(e 3 z)s 

Hence, equation (|3.1I) takes the form 



l + e?k(e?z) s (i + e U(eh) s y 



(3.7) 



R(u) = u zz + u ss + Bi(u) — e 3 u (u 2 — a(e3s,£3z)J = 0, 



(3.8) 

in the region p. 41) . Observe that all terms in the operator B\ have ef as a common 
factor, more precisely we can write 

222 4 22 2 22 

Bi(u) = £3 tt0 (£3 s, £3 z ) Ms + £3 S ai(s3 s,ez z)u z + e 3 sci2(s 3 s,e 3 z)u zz , (3.9) 

for certain smooth functions cijit, 9), j = 0, 1, 2. 

We now consider a further change of variables in equation (|3.8[) with the property 
that it (roughly) replaces at main order the function a by one that has constant 
gradient along T e . Let /3 be as in (| 1 . 36[) . then we define v(x, z) by the relation 

u(s, z) — (3{s^ z)v{x, z), 

(3.10) 

x — (3(s3 z)s. 

Choosing a smaller So, if necessary, we may assume that the coordinates (x, z) are 
also defined for \x\ < 5oE~ 3 , z £ [0,£~3^ e ). We want to express equation (|3.8[) in 
terms of these new coordinates. We compute: 

e^/3 3 v xx , 

\ u z = efl'v + (3v z + ef3'xv X} (3-11) 

eifi'v + 2ei(P') 2 f3- 1 xv x + 2e/3'v z + 2ef3'xv xz + ei /3v zz 

+e ! l(3"xv x + el(l3') 2 p~ 1 x 2 v xx . 
In order to write down the equation, it is also convenient to expand 



/ 1 2 2 2 2 2 4o J 1 .361 „o 2 2 2 4 , 

a(£3 s,e 3 z) = a t (0,e 3 z)e 3 s + a 3 (£3 s, £3 z)£3 s = -(3 e 3 x + a 3 {e 3 s, £ 3 z)e 3 s , 

(3; 12) 

for some bounded function 03 (t, 6). It turns out that u solves (|3.8j) if and only if u, 
defined by (|3.10|) . solves 

S(v) ee £-3^-3^) = Vxx - v (y 2 + x ) + b 2 (v) = 0, (3.13) 

where B2(v) is a differential operator defined by 

B 2 (v) = e^a 3 x 2 v + £"3/3" 3 u zz + e~i /T 3 £i(u), (3.14) 

where the bounded function a,3(t,9) is evaluated at (£^/3 _1 a;, £^ z), and i?i is the 
differential operator in p.7[) where derivatives are expressed in terms of formulas 
(13.111) and s replaced by (3~ 1 x. (Note that B 2 {v) = (3~ 2 v zz + B2(v), where all the 
terms in the operator B2 carry e 3 in front of them) . 
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Our first criterion for choosing an approximate solution v of (13. 13[) is that S(v) 
is small. In view of (|3.13[) . (|3.14[) . it is natural to choose a v that depends only on 
x, and solves the second-order non-autonomous ordinary differential equation: 

v xx - v(v 2 + x) = 0, ieB, (3.15) 

which is known as the Painleve-II equation, a particular case of the second Painlevc 
transcendent [99] . Then, keeping in mind that the inner profile e» (3v((3s) of (|3.1I) 

should match |198| Chpt. 5] with the outer profile J a(e%y) + , as s -> ±oo, it is 
easy to see that the appropriate asymptotic behavior of v should be 

v(x) — \f— x — > 0, x — > — oo; v(x) 0, x oo. (3.16) 

(Observe that the s variable is defined in (— 5oE~i ,5oS~%), but the natural domain 
for this variable is the infinite line). More precisely, recalling (|3.12j) . we get that 



^a e {e'iy)+ = e* [3yf^fs + 0(e\s\i) 



if — 5qE s < s < 0, as e — > 0, while on the other side it holds that a e (e»y) = if 
< s < 5 £-i. 

The following proposition holds: 

Proposition 3.1. Problem \3.1^ )- l3.16\) has a unique solution V, called Hastings- 
McLeod solution. Furthermore, we have 

V x < in M. (3.17) 

The following estimates hold: 

-C\x\~i < V{x)-yf=x< 0, 14 -oo; < V(x) < Cx-ie'^, x>l, (3.18) 
for some constant C > 0, and 

( V x =-l(-x)-i +0(\x\-i), V xx = -±(-x)~i +0{\x\-i), x -> -oo, 

V X V = -± + 0(\x\~ 3 ), V 2 +x = 0(\x\- 2 ), x^-oo, 

\V X \ + \V XX \ <Ce- cx \ x>0. 

(3-19) 

The potential of the associated linearized operator 

M(<p) = Vxx - (3V 2 (x) + x) p (3.20) 

satisfies 

W 2 (x) +x > c> 0, (3.21) 

Proof. We will provide a sketch of proof, underlining the main ideas, and refer the 
interested reader to the original works. 

Existence and uniqueness of a solution V of (|3.15l) - (|3.16p have been shown by 
Hastings and McLeod jl 13J . |116j using a "shooting" technique (they also mention 
an unpublished proof of Conley). Recently, a new proof of existence has been given 
in [18) by the method of upper and lower solutions. Motivated from some problems 
that we mentioned in Section [L2l (see also Remark 13.131 below - ) . we will provide 
in Appendix [B] a new different proof of existence and uniqueness which seems to 
be simpler and more flexible. Another different and rather short proof has been 
given very recently in |2021 Lemma 2.4]. Relation (|3.17[) has also been shown in the 
aforementioned references, and in Appendix [Bl herein. We note that, by looking at 



THE GROUND STATE OF A GROSS-PITAEVSKII ENERGY 



20 



the equation for Z = (— x)~^V, one can show that V has a unique inflection point 
(see (1 13j ) : notice also that V — \J —x is convex whenever nonnegative. 

The asymptotic behavior of V, as x — > ±00, is described in great detail in 
Theorem If .7 of [99] (see also |f 29] ). In particular, it follows that relations (|3.18|) — 
(|3 . 1 9|) hold true. Let us provide some intuition behind these, at first sight, complex 
formulae. The first relation of (|3.1 81) can be formally derived as follows: Assume 
that V — \J—x ~ a(— x)~@ as x — > —00, for some a € K, (3 > 0; plugging this 
ansatz in ()3.15j) . keeping in mind that V + \J —x ~ 2y— x and (the expectation) 
that V xx ~ — j(— ai)~2 as a; — > —00, we readily find that a = — g, /3 = | (see also 
[1761 pg. 160]). A rigorous derivation can be given by simply writing down the 
equation for the function V — \J—x (or better yet for V — \f—x + |(— a;) _ 5) 7 and 
then applying in the resulting identity the following simple lemma, which can be 
proven by a standard barrier argument: 

Lemma (cf. [1041 Lemma 2.1] or |134[ Lemma 3.10]) Suppose that <\> € C 2 , q € C 

-0" + g (x)0 = o(| a; r Q ), 

</> — > 0, and q(x) > c|x|, as x — > —00, /or some constants a, c > 0. T/ien we have 

(j) = 0{\x\~ a ~ 1 ) as x^r -00. 

In passing, we note that the above lemma can also be proven by rewriting the 
equation in terms of the new independent variable £ = — (— x)? and then applying 
a lemma of Herve and Herve [1201 pg. 435], see also |117[ Prop. 3.1]. On the other 
side, because of the second condition in (|3.16|) . the nonlinear term in (13.15|) can be 
(formally) neglected for x — > 00, yielding Airy's equation |32[ pg. 100], namely 

y" = xy, (3.22) 

predicting that 

V{x) - jAi(x) as x -> 00, (3.23) 

for some constant 7 > 0, where Ai(-) denotes the standard Airy's function (see also 
(|3.24[) below). For future reference, we recall from [35] that the two independent 
solutions of Airy's equation can be taken to have the asymptotic behavior 

Ai(a;) ~ i^-^- 1 /4 e -2, 3 / 2 /3 and Bi(x) _ ^^-i/4 e 2- 3/2 /3 j x 00) (3, 24) 

see also the discussion leading to formulae (|B.19I) below. Relation (|3.23[) . which 
clearly implies the validity of the second estimate in fl3.18p . can be established 
rigorously directly from the variation of constants formula 

POO 

V(x) = 7 Ai(x) + 2 / {Ai(x)Bi(i) - Bi(x)Ai(t)}V 3 (t)dt, (3.25) 

J X 

see also [16j Lemma B.l]. It is worthwhile to note that the exact value of 7 was 
determined to be 

7 = 1 (3.26) 

in [113] . using an integral equation derived by inverse scattering techniques of 
Ablowitz and Segur. An implication of this beautiful formula, for the physics of 
Bose-Einstein condensation, has been pointed out in [161] . The estimates in (|3 . 19[) 
can be established in a similar manner. For example, the asymptotic behavior of V x 
as x — > —00 follows by writing down the equation for the function V x 4- |(— ic)~2 
and then applying the above lemma. 
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Finally, the lower bound (|3 . 2 1 1) lias been proven recently in 11)1. Lemma 2.2]; 
despite of its simple appearance, its proof takes almost three pages! Actually, 
we were surprised to find that such an apparently simple relation proves to be so 
recalcitrant. We emphasize that, for the purposes of the current paper, it can be 
bypassed (something that we did not do) at the expense of some extra work (see 
Remark |3 . 71 below - ) . On the other side, it is quite painless to show that the operator 
M. still satisfies the maximum principle without knowledge of the latter lemma. In 
other words: 

"whenever w € C 2 (R), M.{w) < 0, and liminf w(x) > 0, then necessarily w > 0". 

\x\— >oo 

(3.27) 

For this, recall that the typical way towards establishing the maximum principle for 
an elliptic operator L in an open set f2 is to first show that 

"there exists a positive upper-solution (j> of L(cf>) — in Q", (3.28) 

see [51] for more details. If this holds, then adding one of various additional as- 
sumptions on cj> (the simplest one being <f> > c with c a positive constant), it does 
guarantee the maximum principle to hold, see |35j . Indeed, in bounded domains, 
(13.281) is a necessary and almost sufficient condition for the maximum principle to 
hold, see Corollary 2.1 of [35]. However, in unbounded domains one has to be more 
careful to deal with infinity. Having this in mind, firstly note that 

M(-V x ) = —V < 0; M(V) = -2V 3 < 0, (3.29) 

and recall that -V x > 0, V > in R. So, condition (J3T28]) . with L = -M, is 
satisfied by either one of —V x or V. Nevertheless, observe that one faces a difficulty 
when proceeding as in |35j , namely applying the standard maximum principle in the 
equation satisfied by ^y- or y (where w is as in (j3.27jl ). This is because —V x and V 
vanish as x — > ±oo and x — > oo respectively, making the signs of liminf x ^.±00 -^rr 
and lim infx-yoo y unclear (in the case where one of liminf ^^±00 w equals zero). 
(This subtle point seems to have been overlooked in 202 ). Instead, we consider 
the function 

w-eV x 
<Pe = y , e > 0, 

noting that, thanks to (|3.22p - (|3.24[) (see also (|4.25|) below), we have lim^^-oo ip e > 
and lim x _ ) . 00 ip c = 00. Using (|3.27j) . (|3.29|) . and applying the standard maximum 
principle in the equation satisfied by ip e , we obtain that ip e > in R, for every 
e > 0. Then, letting e — > we infer that w > in R which is the desired assertion 
in (13371) . 

The sketch of the proof of the proposition is complete. □ 



We can now define the inner solution of problem (13.11) . in the neighborhood of 
T £ described by {|x| < S £~~^, z e [0,£~3£ e )}, as 

u m {y)=eif3(eh)V(x), (3.30) 

(recall (|3~TU1) '). 

Let L > be a large constant to be determined independently of small e > 0. 
We consider the neighborhood of r e described by 

U e = {-2L < x < 5 e-^, ze[0,s-i£ s )}. (3.31) 

The following proposition contains the main estimate regarding Uj n . 



THE GROUND STATE OF A GROSS-PITAEVSKII ENERGY 



22 



Proposition 3.2. If e is sufficiently small, the inner approximation Ui n , defined 
in L3. 30\) . satisfies 

\\Au m - £-iu in (uj n - a{eiy)^ \\ L -{v.) < Ce. (3.32) 
Proof. From the calculation leading to (|3.13[) . and (|3.14p . we find that 
Au in -e-iu in (u 2 in -a(eiy))\ = \R(u m )\ = ei[3 3 \S(V)\ 



■ Ce 



(3.33) 



(x 2 +es)V+ \V X \ + \x\\V xx \ 



pointwise in U e . The desired estimate (|3.32[) now follows at once from the above 
relation, via (f3TT8|) . ff37T9]t and ([331]) . 

The proof of the proposition is complete. □ 

Remark 3.2. From the geometric singular perturbation viewpoint [190 \, the Hastings- 
McLeod solution corresponds to a special trajectory of the (de-singularized) blown- 
up system that connects two equilibria on a sphere. The fact that this connec- 
tion is transverse, which allows for a perturbation argument, follows from the non- 
degeneracy of the Hastings- McLeod solution (as defined in Remark \3.7\ below). In 
this regard, see also Remark \3.5\ below. 



Remark 3.3. In the context of singular perturbation problems, the Hastings-McLeod 
solution first served as a basis for an inner solution in plasma physics (see [76] ). 
Since then, it has been (formally) used to describe layered structures in problems 
involving crystalline interphase boundaries [TH 11901 1191] . patterns of convection 
in rectangular platform containers [75], self- similar parabolic optical solitary waves 
[41j , and the Navier-Stokes and continuity equations for axisymmetric flow [157j . 

Remark 3.4. In well known singularly perturbed elliptic problems, such as the 
focusing NLS (recall the discussion in Subsection M .2.2\) or the spatially inhomoge- 
neous Allen-Cahn equation (see for instance 86J, [151] . [201] ), the corresponding 
inner profile is determined by special solutions of autonomous second order elliptic 
equations (posed in less or equal dimensions). In the former case the corresponding 
equation is 

Au-u + \u\ q ^u = 0, 

while in the latter it is equation l[1.4-9\ ). In contrast, in the problem at hand the 
corresponding equation &3.15\) is non-autonomous, as was the case in [37|, [TT] . 
[73j , [134] , and [135] . The interested reader can verify that similar situations also 
occur in the singularly perturbed Fisher's equation [1 181 Chpt. 10, Exc.3] (see also 
[165] ). and in the spatially inhomogeneous Allen-Cahn equation that we mentioned 
previously, treated in [86] , [151] , if the spatial inhomogeneity is not strictly positive 
but vanishes at certain points ( or submanifolds ) of the domain (have in mind the 
first part of Subsection M.2.3} ). 

3.2. Set up away from the curve. In this subsection, adapting an idea of |134] 
(see also [135J), we suitably perturb, in T> e , the outer approximation 



Uout 



\/a(eiy)+ (3.34) 



to an improved outer approximation u out , which is closer to the inner approximation 
Ui n near the curve T e , for small e > 0. We emphasize that this is a key step 
in our construction of an approximate solution for ([3.1jl . Our underlying plan, 



THE GROUND STATE OF A GROSS-PITAEVSKII ENERGY 



23 



carried out in the following subsection, is to smoothly interpolate between u out 
and Uj, n , and Ui n and zero (in M. 2 \T> £ ), near the curve r e , in order to obtain an 
approximate solution u ap of (|3.ip that is valid in all of M 2 . Interpolating directly 
between u ou t, as defined in (|3.34p . and Ui n in T> £ is not standard due to the following 
obstructions: The inner approximation Ui n leaves a remainder in fl3.1[) that grows 
with respect to the (negative) distance from r e , as can be seen from (|3.18[) , (|3 .33|) : V 
converges algebraically slowly to \J— x asu — oo, recall (|3.18|) . (For more details 
on this subtle point we refer to Remark 13.81 below) . In the one-dimensional case, an 
elegant solution to this can be given by geometric singular perturbation theory, see 
[190] and Remark [X5] below. In contrast, thanks to the (super ) exponentially fast 
convergence of V to zero, as x — > oo, interpolating directly between m n and zero in 
K 2 \2? £ is rather standard. (This last situation occurs in the construction of spike- 
layered solutions of the focusing (ll.22[) + . see [5D], HE], [SE]> an d transition-layered 
solutions of Allcn-Cahn type equations, see [75], [55]). 

Let 5 < ioo(i+max/3) a nxe d number. We consider a smooth cutoff function 




n 5 {t)={ ( :-!.:-!■-,) 

Denote as well 

Xs(x) = n s (e^x), (3.36) 
where x is the coordinate in (|3.10[) . 

We define our outer approximation in the domain T> e \{— L < f3s < 0, z G 
[0,£"i4)} to be 

1 

u out {y) = [a(eiy) + s%(3 2 Xs (p S ) [/3s + V 2 (f3s)] }' . (3.37) 

Note that, thanks to (|1.33j) . (|3.19|) . if L is sufficiently large then u out is well defined 
for small e. 

The following proposition makes delicate use of estimates (|3.18j) . ()3.19|) . and 
contains the main properties of u out . In some sense, it is the "heart" of the present 
paper. 

Proposition 3.3. We have 

Au out - £-iu out (u 2 out - a(s%y)j = 0(e|s| - ') (3.38) 

uniformly in {— Se~i < j3s < —L, z € [0, s~i£ e )}, and 

Au out - s'iuout (u 2 out - a(e*y) \ = 0(e*) (3.39) 

uniformly in f> £ \{~ 5e~i < (3s < 0, z € [0, e~il e )}, as e — > 0. Moreover, if e > 
is sufficiently small, we have 

\u ou t - u m \ < Ce\s\^ , \V(u out - u in )\ < Ce\s\^ , \A(u out - u in )\ < Ce\s\~^ 

(3.40) 

in {— 5e~3 < (3s < —L, z e [0,e _ 3 C-s)}, and 



Uout - \Ja(ely) 
m {-28e-§ <(3s<~L, ze [0,£-§4)}- 



<Ce*|s|-* (3.41) 
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Proof. In {-de a < (3s < —L, z 6 [0, e ^£ e )},we have xs — 1 and we can compute 
that 



-e lu out (ul ut - a(ely) 
via Ol . (OoD 



via (l3~T5)) 

via (l3~18l) . (l3~19l) 

uniformly, as e — > 0. Moreover, in the same region, we have 



= — e 



-/3 2 a{eiy) + si^s + s^p 2 V 2 ' [(3s + V 2 {(3s)} 

-(3 2 \o(eh 2 ) + ei(3 2 V 2 } ' [f3s + V 2 (f3s)] 

-e^ 3 V((3s) \l + elO{s 2 V~ 2 )} 2 [j3s + V 2 {Ps)] 
-e^(3 3 V{(3s) [l + el0(|s|)l [(3s + V 2 ((3s)] 
^ 3 V((3s) [(3s + V 2 ((3s)] +0(e\s\-i), 



(3.42) 



(Uout)i 



a t (eh,sh) + p 3 + 2(3 3 V x {(3s)V{(3s) a{s%y) + e%(3 3 s + ei(3 2 V 2 ((3s) 

(3.43) 

£§a t4 (£§s,£iz) + 2/3 4 (V xx V + V 2 ) a{siy) + ei(3 3 s + ei (3 2 V 2 



Eg 

2 



g5 
4 



a t (eis,eiz) + (3 3 + 2(3 3 V X V a{e^ y) + (3 3 s + (3 2 V 2 

(3.44) 



(u ou t)z 



and 

(Uout)z 



e^ae + ^(3'(3 2 s + 2e^ (3' (3V 2 + 2e^(3'p 2 sV x V 



a + ei(3 3 s + ez(3 2 V 2 
(3.45) 



:ta ee + 6e 2 (P') 2 (3s + 3e 2 P"(3 2 s + 2e 2 V 2 (p"(3 + ((3 1 ) 2 ) + 8e 2 (P') 2 (3sV x V 



+2e 2 (3"(3 2 sV x V + 2e 2 ((3') 2 (3 2 s 2 (V 2 + V XX V)} 



a + e^(3 3 s + ei(3 2 V 2 



-\ esag + 3s3(3'(3 2 s + 2e3(3'l3{V 2 + (3sV x V) a + (3 3 s + (3 2 V 2 

(3.46) 

where the functions a, ag aee are evaluated at £$y, V at (3s, and (3 at e^z. Esti- 
mating as in p. 421) . we conclude that in the region described by {— 5e~3 < (3s < 
—L, z e [0,e~3£ £ )} we have: 

(u out ) s = e^ 2 V x ((3s) + 0{e\s\i), (u ut)ss = (3 3 V xx ((3s) + 0(e\s\-i), 
{uout)z = C(e|s|5), {uout)zz = C(e|s|~"5 ), 

(3.47) 

uniformly as e — > 0. We point out that, when estimating derivatives in z, we also 
made use of the bound 

\ag(es s,£3 z)\ < Cez |s|, 
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which follows directly from (| 1 . 33|) . In order to elucidate the various cancelations of 
powers of \s\ involved, let us carefully present the proof of one of the estimates in 

\(u ou t) zz \ < C(ei +e 2 \s\ + £ 2 | S ||s|-5| s |l +s 2 .s 2 (\ S \- 1 + M"*M*)) e^M"* 
+C(e 



^£3 \S\ + £3 \S\\S\ 2 \S\ 2 j £ \S\ 2 

< C ^ +s 2 \s\^s-i\s\-^ +Cel\s\ 2 e~ 

< Ce\s\-i +C£§ |s|l < Ce|s|-3 +C£l|s||s|-3 < Ce|s|-4. 

In {-2&ri < /3s < -5e _ §, z € [0,£~f4)}, making again use of (f3TT8|) . (f37T9|) . 
we can show that 

w 2 uf - a(eiy) = 0(s 2 ), (« out ) a = 0(ei), {u out ) ss = O(e^), 

(3.48) 

(Sout)z = 0(£3), (u „t) 2Z = 0(es), 

uniformly as e — > 0. Estimate (|3.38p now follows from (|3.6[) . (|3.9p . the fact that V 
solves ([37T5]) . ([3TT9"]) . (|3~22]) . and (j3T471) ; estimate (|3~351l follows at once from 
(|3.9p , (|3.37[) , and (I3.48P . The first relation in (|3.40[) can be shown by working as in 
(1552)) , The remaining estimates in ([3~I0"j) follow directly from (l3TTg|) , (133^]) , ([330)1 , 
(j3T47j) . and the fact that 

u yi + U L = u * + (1 + ehs) 2 u 2 z , (3.49) 

for y = (2/1,2/2) € {(s,.z) | |s| < <5o£~^, z € [ti,s~i l E )}, and any smooth function 
u defined in this region (recall Q3.5P ). Finally, estimate Q3.4ip follows readily from 
(PUD, and dSHSJ). 

The proof of the proposition is complete. □ 

3.3. The approximate solution u ap . In this subsection we will construct a smooth 
approximate solution for problem (|3.1[) that is valid in all of M 2 . This will be 
achieved by smoothly interpolating between u Q ut and Ui n in T> £ , and between Ui n 
and zero in M. 2 \T> £ . 

To this end, we need to introduce one more smooth cutoff function: 

( 0, x > —L, 
Pl(x) = \ (3.50) 
[ 1, x < -2L. 

We can now define our approximate solution for (I3.ip as 

' u out in T> S \{-2L < x < 0}, 



u a P — \ u in + p L (x)(u out ~ u ln ) in {-2L < x < 5e I}, (3.51) 

k Xws(x)u ln everywhere else, 

where Ui n , \5i u out were defined in (|3.30p , (|3. 361) . and (13.371) respectively. 
The following proposition contains the main estimates concerning u ap . 
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Proposition 3.4. The approximate solution u ap satisfies 



ap 

Au ap - e^u ap (u 2 ap - a(e*»)) = O (e(\s\ + 1)-*) (3.52) 
uniformly in {— 5e~i < (3s < 0, z € [0,e _ 3£ e )} ; 

Au ap - e"3 Uap ^ _ a ( e f = 0( £ |) (3.53) 

uniformly in T> e \{—5e~i < /3s < 0, z £ [Oje - ^,^)}, 

Au ap - e-iu ap (u 2 ap - a(eiy)) = 0{ee- c ^) (3.54) 

uniformly in {0 < (3s < 28e~i , z € [0,£~^i? £ )} 7 as e — > 0, and 

u ap = everywhere else. (3.55) 
If e > is sufficiently small, we have 

{ces(l + |a;|), «/ l^l^fe - ^ 
(3.56) 
c + c\e^y\ p , otherwise. 

Proof. We will first consider relations p.52j) - (l3.55[) . In view of estimates (|3.32[) , 
(|3.33p . (I3.38p . p.39p . and recalling the super-exponential decay of V as x — > oo, it 
just remains to show the validity of f|3 . 52[) in the interpolating region described by 
{-2L < (3s <-L,zE [0,£-i£ e )}. There, we have 

Au ap - e~iu ap (u 2 ap - a(etj/)J = Au in - s~^u in {uj n - a(s%y)\ + (Ap L )(u out - u in ) 

+2Vp L V(u out - U in ) + PhA(u out - Uin) 
-£~iu in [ip 2 L {u out - U in ) 2 + 2u in p L (u out - U in )] 

s~^pL{ucmt ~ u in ) u 2 n - a(ei y) + p 2 L (u out - u ln ) 2 

and the desired estimate follows via p.32[) . p.40[) . noting that uf n , a(eiy) are of 
order £3 in this region. 

The proof of lower bound (|3.56|) proceeds as follows: In the neighborhood de- 
scribed by {\x\ < 2L} of the curve T e (with the obvious notation), by virtue of 
(f03) . dmDJl, (|33T|) . we find that 

3u 2 ap -a(eiy) = iei(3 2 V 2 {x) + 0{e$) + e%(3 2 x + 0{fJ)x 2 

(3.57) 

= eip 2 (W 2 (x)+x) +G(e§) 

uniformly as e — > 0. In T> e \{— 2L < x < 0}, we have u ap = u out and, by ()3.19[) . 
p.37|) . we infer that 

3ul p -a(eiy) = 2a{eiy) + 3ei/3 2 X5 (/3s) (V 2 (p S ) + /3s) 

(3.58) 

= 2a(eiy) + eh S ((3s)0(\s\- 2 ) 

uniformly as e — > 0. In points outside of the domain T> e U {0 < x < 2L}, we plainly 
note that 

3u 2 ap -a(eiy)>-a(eh)- (3.59) 
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The desired lower bound (13.561) now follows readily from the above three relations, 
via ()1.4|) . (|1.33[) . and (|3.21[) . increasing L > if necessary. 

The proof of the proposition is complete. □ 

Remark 3.5. From the geometric singular perturbation viewpoint, recall Remark 
YS.'A matching is accomplished by employing a useful lemma on the flow past a "cor- 
ner equilibrium" (see [189U - Manifolds of corner equilibria arise in blown-up geo- 
metric singular perturbation problems precisely where the inner and outer solutions 
must be matched. When such equilibria are normally hyperbolic, as in the one- 
dimensional case of the problem at hand (see |190j ). this lemma plays the same role 
in tracking the flow past them that the Exchange Lemma [131] plays at certain other 
manifolds of equilibria. 

Remark 3.6. Our construction ofu ap should also be applicable to the homogenized 
Gross-Pitaevskii equations considered in [162( Sec. 7] . 

3.4. Mapping properties of the linearized operator. In this subsection we 
will invert the linearized operator 



C(<p) =A<p- e-S (3u 2 ap - a(eiy)) f 



(3.60) 



in carefully chosen weighted spaces. The use of weighted spaces is a powerful tech- 
nique in elliptic singular perturbation problems, and in many problems arising from 
geometry, see |171j . To the best of our knowledge, they are used here for the 
first time in singular perturbation problems involving corner layers. Actually, the 
weighted spaces that we will use are a variant of those considered in [179j , and are 
motivated from relations (|3.38p . (|3.56[) . keeping in mind that we ultimately wish to 
find a true solution of (|3.ip near u ap via a perturbation argument. 
Consider a smooth non-increasing function g such that 



1. 



9(s) = < 



max (3 \ 2 

L J \ max /3 



s > 0, 



2t 



min /3 



£ 3 < S < 



(3.61) 



i max P \ 



L 2£ 



s < 



36 
min j3 



£ 3, 



and 



< -g' < CL~ X , \g"\ < CL,- 2 , s € 



L 



max f3 



.0 



o < -g' < CL~2£-3, \g"\<CL-^£t, se 



36 
min/3 



£ 3 



26 
min/3 



£ 3 



(3.62) 



(3.63) 



where the constant C is independent of small e and large L. Recalling (13.56[) . it is 
easy to check that we can fix an Lq > such that 

_2 

< Ce^ + CL- 2 < E —^- min [zu\ p - a{e*y)) , s G R, if L > L , 

(3.64) 

provided e is sufficiently small (C in the above relation is independent of e, L). 
Relations (13.571) — (|3 . 50[) imply that, for small e, we have 





+ 2 


I 


9 




9 



9*(y)<Ce 3 (3ul p -a(£3y)) in 



(3.65) 
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(here g is viewed as a smooth function of y, which close to T £ , in coordinates (s, z), 
is given by (|3.61|) . and otherwise equals the constants in (|3.6ip ). 
For p € L°°(R 2 ), we define the following weighted norms: 

IMIf = ll3^lli~(R 2 ) and \\<p\\i = \\g*<p\\ L °°(m 2 )- (3-66) 

(We utilized this notation because g behaves qualitatively like (— s)i for s < 0). 
We also consider the Banach space 

X = {p : \\p\\ x = || e l £§? '^|| L oo (R2) }<oo. (3.67) 
The following proposition will be used essentially in the sequel. 

Proposition 3.5. If e is sufficiently small, given f 6 X n C a (R 2 ), < a < 1, 
there exists a unique p 6 X l~l C 2+Q (K 2 ) such that 

C(p) = f, (3.68) 

where the linear operator C was defined in \3. 60\) . Furthermore, we have 

IMk<C||/|U, (3.69) 

and 

IMI|<C||/||i, (3.70) 
for some constant C independent of f,e. 

Proof. The first assertion of the proposition, including estimate (|3.69p . follows in a 
standard way: It follows from (|3.56[) , the maximum principle, and elliptic regularity 
theory [TU2], that there exists a solution ip e C" 2+Q (IR 2 ) of ([335)1 such that 

IMIl~(r 2 ) < C||/|| L ~( R 2), (3.71) 

for some constant C independent of /, e. This is easy to prove, though it is difficult to 
find a good reference. (For example, one can first solve equation (|3.68[) in a ball Br 
with Dirichlet boundary conditions to obtain a solution pn such that ||¥?.r||.l° o (.Br) < 
C||/IU»°(B K )> for some C independent of f,e,R, extend pr in L°°(R 2 ) to be zero 
outside of Br, and prove that p^ , for some B4 — > 00, converge uniformly in compact 
sets of K 2 to a solution ip of (|3.68[) that satisfies (|3.71[Q . Then, a standard barrier 
argument, using as barrier the function 

MV) = + IMU^e^-le^l), \ y] > J_ j (3 . 72) 

where R > is chosen large, r > arbitrary, yields that \p{y)\ < <p~r(y)> \y\ > 
i?e~ 2 / 3 , provided e is small enough so that (|3.56[) holds. Lastly, letting r — > 0, and 
recalling (|3.71|) . we conclude that estimate (|3.69j) holds true (see also [8J3 Lemma 
7.3]). 
Let 

tp = gp, (3.73) 
then, thanks to (|3.6p . p.9p . it is easy to see that, with the obvious notation, we 
have 

A^-2^ s -£-f (Zul p -a{eiy))ip- 9 jtp + 2(^j tp-e^a Q 9 -^ = gf, (3.74) 

y £ M 2 . Since -0 — >■ as |y| — 00, recalling (|3.56[) , (|3.64p . p.65p . we can apply the 
maximum principle to show that, for small e, 

\Tp(y)\ <C||^/|U» (R2)) y eR 2 , 
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for some constant C independent of e, f. We also used the fact that ips = whenever 
Vyip — 0, which follows immediately from relation (|3.49|) . The desired bound (|3.70| 
now follows at once from ()3.66p . (|3.73[) . and the above relation. 

The proof of the proposition is complete. □ 

Remark 3.7. In the above proof, we made essential use of lower bound US. 56]) whose 
proof, we recall, relied crucially on lower bound 13.21]) which was established recently 
in [104j . However, as we have remarked in the proof of Proposition \3.1[ its proof is 
rather involved and technical. In AvvendixVfi[ we will provide a more natural and 
flexible proof of Proposition Iff. 51 without assuming knowledge of i3.21]) . Instead, 
we will make use of the asymptotic stability of the Hasting s-McLeod solution V , 
namely the fact that the principal eigenvalue of the operator —M., defined in (3.20]) . 
is strictly positive. This follows immediately by testing the corresponding eigenvalue 
problem by V x < 0, see |191j and Proposition \B.1\ herein. In other words, we will 
rely on the lower bound: 

/CO 
(fdx V^C M (R), (3.75) 
-oo 

for some constant c > 0, which clearly is much "softer" than \3.21]) . We point 
out that the validity of (3. 21]) was not needed (nor known) in references [190] , [191] 
which dealt with related one- dimensional problems. Actually, in one-dimensional 
or radially symmetric cases, for Proposition Iff. 51 to hold, it suffices to know that 
zero is not in the kernel of M. . However, this may not be true in general higher 
dimensional problems due to a possible resonance phenomenon (see [78], [135] ). 

3.5. Existence of a solution. Here we will use the contraction mapping principle 
in order to capture a genuine solution u e of (|3.1|) close to the approximate solution 
u ap , provided e is sufficiently small. 



Proposition 3.6. If e is sufficiently small, then there exists a solution u e of 13.1]) 
such that 

IK - u ap \\x < Ce, (3.76) 

K-u«p||j <Ce, (3.77) 

and 

\u s (y) - u ap (y)\ <Cee~ c \ y e {0 < /3s < 2fc~t, z G [0,£-§4)}, (3.78) 
where the norms involved were defined in i3. 66]) and |ff.6'7| ). 
Proof. We seek a true solution of problem (|3.1[) in the form 

u — u ap + (p. 

In order for u to satisfy the equation in (|3.1|) . we readily find that the correction ip 
has to solve 

C(cp) = E + N(<p), (3.79) 
where the linear operator C was defined in (|3.60[) . 

E = -Au ap + £~iu ap (u 2 ap - a(£%yf) , (3.80) 

and 

N(ip) = 3e-iu ap <p 2 + e-i^ 3 . (3.81) 
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Given M > to be determined (independently of e), we consider the closed 
bounded subset of X denned by 

y = {<peX : \\<p\\ x <Me, |M||<Me}. (3.82) 

We will prove that, if M is chosen sufficiently large, the operator P, defined from 
X n C 2+Q (R 2 ) into X n C 2+a (R 2 ) by 

f( (P ) = c- 1 (e + n( v >)), 

maps y into itself, and is a contraction with respect to the A"-norm, provided e is 
sufficiently small. Note that £ _1 is well defined by virtue of Proposition 13.51 In 
view of the estimates of Proposition ^. 4[ (|3.6ip . (I3.66|) . and (|3.67p . for small e > 0, 
we have 

\\E\\ X < Ce and ||E||i < Ce. (3.83) 

Furthermore, there exists a constant C > such that, for all tpi,ipz,ip G X, the 
following relations hold pointwise: 

\N(<p!) - Nfa)] < Ce~^((pl + ^)|<^i - pa| + Ce~i \u ap \{\<pi\ + " M 

\g*N(<p)\ < Ce-^g^l 3 + Ce-^\u ap \\gip\ 2 , 

(3.84) 

for every y G R 2 (recall that g > 1). If y> G 3^, by Proposition 13.51 (|3.83l) . and 
(|3.84|) . we obtain that 

\\P(ip)\\ x < C\\E\\ x + C\\N(tp)\\ x 

< Ce + Cs-§M x +Ce-§y\\ x 

< Ce + CM 3 ei + CM 2 e^ , 

and 

\\p( v )\\ 3 _ < cp^ + qiJvfoOiij 

< Ce + Ce-I||^|| 3 3 +Ce-i||^|| 2 

2 2 

< Cs + CM 3 e% +CM 2 ei, 

where C is independent of <p,M and small e. We conclude that, if M is chosen 
sufficiently large, the operator P maps y into itself, provided e is sufficiently small. 
We have to prove that P is a contraction from y into itself with respect to the 
A"-norm. Let tpx^tpz € y. As before, we have 

||Pfo>i)-Pfo> a )IU < C||A%>i)- AWIU 
< Cei||</3i - </3 2 |k- 

Hence, for small £ > 0, the operator P : y ^> y becomes a contraction with respect 
to the A"- norm. So, recalling that y is closed in X, it has a unique fixed point 
<p* G y, thanks to the contraction mapping theorem (see for instance |109j ). It is 
clear that the function u e = u ap + (p* satisfies the elliptic equation in (|3.1|1 . and 
estimates (j37T6]l - (prF?l) . 
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Next, we show that u £ is positive, and consequently solves problem (|3.1[) . In 
the neighborhood described by {|a;| < 2L} of the curve T e , recalling (|3.30p , (j3.40|) . 
(|3.5ip . and that ip* e y, we have 

u £ = Mm + O(e) = e^[3V(x) + 0(e) > ce^ 

uniformly as e — > . In the domain T> £ \{— 2L < x < 0}, thanks to (|3.19j) . (|3.37|1 . 
we have 

2 2 1 3 11.331 1 

u E = Uout + 0(e) = a(e3y) + 0{e3L- 2 ) + 0(e) > ce^ 

uniformly as e — > (having increased L if necessary) . It remains to consider points 
outside of the domain T> e U {0 < x < 2L}, where u solves an equation of the form 

Au — p(y)u = 0, where p > c, (3.85) 

(recall (|1.4p . (|1.33[1 ). The positivity of u in this region follows directly from the 
maximum principle, and the fact that we have already shown that u > ce^ on the 
boundary described by the closed curve {x = 2L}. 

It remains to establish the validity of (|3.76l) - (|3.78p . Since u e — u ap = cp* e y, we 
see that (|3.76j) . (|3.77[) hold. Finally, we will show (|3.78|l by suitably modifying the 
proof of Lemma 2 in [38] (see also [95l Pg- 230] and [1251 Lem. 2.2]). From (|3T54|) . 
(|3"35f . (j3~56l) . (j3~76l) . and (f3~79l) . we find that ^* satisfies 

Atp* - P{y)ip* = O (ee~^ s ) , where P(y) > c, (3.86) 

uniformly in {0 < /3s < 4fe~i , z e {0,s~%£ e )}, as e — > 0. The reason for choosing, 
in the righthand side, a decay rate strictly less than yfc is to facilitate our next 
argument. Let 



y € {0 < /3s < 4fe~3 , z £ [0, e~z£ e )}, where the value of the large constant M > 
will soon be fixed independently of small e. By virtue of (|3.76p . (|3.86l) . we can 
choose a large M > such that 

-A(^-</3*) + P(y)(£-(^) >0, ye {0 < ps < 45e~*, ze [Q.e-f^)}, 

fi-tp t ,>0 on {s = 0} U {/3s = 4fc"i}, 

if e is sufficiently small. Now, by the second estimate in (I3.86p . and the maximum 
principle, we deduce that 

<p- <p* > if < (3s < 4fe"i. 

In turn, the above estimate readily implies the validity of (|3. 781) . 

The proof of the proposition is complete. □ 

The following estimates hold: 
Corrolarry 3.1. The solution u e of \3.1\l , constructed in Provosition \3.b\ satisfies 
u e = ei/3(siz)V (f3(eh)s^ + 0{e\s\i + e) (3.87) 
uniformly in {— Se~i < (3s < 0, z G [0, e~i£ e )}, 

u e = ei(3{eiz)V {p{e^z)s^ + 0(e e - cs ) (3.88) 
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uniformly in {0 < j3s < 2fe~i , z G {0,s~%l e )}, 

u E = yJa(e$y) + 0(e*\ a \-%) (3.89) 
uniformly in {— fe~t < /3s < — 2L, z G [0,e~^£ £ )}, 

M e = \Ja{eiy) + 0(e 2 ) (3.90) 
uniformly in T> e \{—8e~i < /3s < 0, z G [0,£~^ e )} 7 and 

< u £ < Ce^ exp{-cdist(y, V e )} (3.91) 

in R 2 \V £ , as e ->■ 0. 

Proo/. Estimates (|3"lT7j) . (g^g| follow readily from ([3~50"j) . the first relation in (1TT4T))) , 

(133T1) . (j3T7S)) . and (j3~78")) . From ([3~4T]) . (j33T|) . (|3~FT|) . (j3Tf|) . for small e > 0, we 
have 



\Ja{ely) 



in {— 5e % < f3s < —2L, z £ [0,e ^£ e )}, and estimate f|3 . 89[) follows immediately. 
Estimate (j3T9U|) follows at once from ([337]) . (|33T|) . (pT6TTl . and ([3J7|) . Finally, 
estimate (13.911) follows readily from (13.85)) , (13.881) , arguing as we did for the proof of 
(13.781) (but here we need to cover R 2 by a finite number of disjoint annular domains, 
surrounding T> , and the exterior of a large ball), see also [95] pg. 230]. 

The proof of the corollary is complete. □ 

Remark 3.8. Notice that the nonlinear terms in \3. 79\ ) are of cubic-like order. 
Indeed, as in the proof of Proposition \3. 51 using the maximum principle, 13.56V . 
and the easily derived bound 

\u ap \ < Ce~i (3u 2 ap - o(e*y)) , y G K 2 , (3.92) 

we can show that, if e is small, the unique solution ip £ X n C 2+a (R 2 ) of 

C(<p) = u ap f, f G X n C a (R 2 ), < a < 1, (3.93) 

satisfies 

IMU < Ce*\\f\\ x . (3.94) 
Therefore, in order to successfully apply the contraction mapping principle, as we 
did in the proof of Provosition [3~6\ it is enough to construct an approximate solution 
v ap supported in a ball of radius 0{e~~s) such that 

£ = — Av ap + e~3v ap fv^p — a(e^y)J = OQ hxe^e^), uniformly in R 2 , as e -» 0, 

(3.95) 

for some constant 7 > 0, and \3.5b}) . \3. 92]) remain true, with v ap in place of 
u a p, for sufficiently small e (the logarithmic term in \3. 95\) is used for convenience 
purposes only and has nothing to do with that appearing in HI. 25)) ). This was the 
main strategy followed in [191] for a related one- dimensional problem. Actually, one 
can plainly define an approximate solution for \3.1\) as 

u out in V e \{-2M E < x < 0}, 

Ui n + PM E (x)(u out - Um) in {-2M e < x < 8e~i}, (3.96) 
Xs{x)ui n everywhere else, 
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where M F is such that L < M F < j^e 3, and Ui n , u ou t, pu as in S3.30\) , \3. 34\ ), 



Ii3.50\) respectively. Working as in Proposition \3.3l we can verify that 

\uout - u in \ + \s\ \V(u out - u in )\ + s 2 \A(u out - u in )\ < C (e\s\% +e^\s\~^ 



in the region described by {— 25e 3 < j3s < —L}. Then, as in Proposition \3.4\ we 
can show that £ , defined in \3.95}) , satisfies 



£ 3 , 



£3| S 2, 



\£\<Cl 



e s 2 + £3 s 2 . 



M £ 2 e + M £ 2 e§. 



m V E \{-he-i < f3s < 0}, 
-fe-§ <f3s< -2M £ , 
-2M E </3s< -M e , 

-M e < [3s < 0, 
everywhere else. 



Consequently, 



achieve bound \3. 95\) by plainly choosing M e — |lne|. Fur- 



the 



rmore, the approximation v ap is sufficiently close to u ap so that the estimates 
56)) and US. 9ty) remain true with v ap in place of u ap . However, the correspond- 



ing estimates for the solution of \3.1\) . obtained using this approximation, are far 
from optimal. One can actually check that the above argument works because the 
exponent 5/2 in \3.18]) is strictly larger than one. It is worthwhile to mention that 
the geometric singular perturbation approach in [190 required merely VS. lb}) . In 

[202] , [203] . for a closely related problem to (jl.22[) . the authors made the choice 

M e = j^£~3 (according to our notation) while at the same time not using any 



convergence rate of V(x) to 
us. 



-x as x — > —00, something which is not yet clear to 



3.6. Proof of the main theorem. We are now ready for the 

PROOF OF THEOREM [TTT1 It follows from the definition of a e from dL32l) that 



u e (y) 



ye 



(3.97) 



where u £ is the solution of p.l[) as in Proposition ^. 61 is also a solution of problem 
(|1.8[) besides the minimizer r\ e of G £ in H. On the other hand, we know from 
Theorem 2.1 in |125] that (given A E ) problem (|1.8|) has a unique solution (see also 
Remark 13.121 below) . Therefore, we conclude that u e = r\ e . Estimates (|1.38|) - (|1.42|l 
for rj £ follow readily from the corresponding estimates p.87|) - (|3.91j) for u £ . Relation 
(fOtj) follows easily from (j33Sl) and (IBTTo) . 

Next, we will derive estimate (|1 .45[) by building on estimates (|1.25|) . (I3.87P 
(|3.9ip . and using that ||??e||l 2 (k 2 ) = 1- (F° r a self-contained proof of (|1.8p we refer 
to Remark |3 . 141 below) . We consider the following annular regions of the plane: 



{-Se-i < f3s< -2L, 0<z< e~^£ £ }, 



So = {-2L <(3s<0, 0<z< e-h £ }, 
S+ = {0 < f3s < 5e-%, 0<z< E~i£ E }. 
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It follows from (|3.51| , and (|3.77J) . that the solution of the stretched problem (|3 . 1 1) 
satisfies 

u £ = u out + 0(e\s\~z), uniformly in S-, as e — >• 0. 
Furthermore, from (| 1 .33[) . f|3 . 19[) . and (|3.37p , if e is small, we have 

csi\s\i <u ut<Ce^\s\i in S-. 
So, from (11.91) . (|1.32[) . (|3.3T[) . and the above two relations, we find that 
u\ - A{eiy) = \ £ - A + eip 2 [/3s + V 2 (/3 S )} + <D{e? Isp 1 ), 
uniformly in SL, as e — > 0. Thus, via the identity 

f(y)dy= / / a f(eh,eh)(l + sh £ (eh)s)dsdz V/ e C(S_), 

'S_ JO J-5s _ 3£-i v ' 

we obtain that 

J s _ (ul - A( £ ly)) dy = (A £ - A ) / s _ Idy + ft (3 £ (8)d8 [V 2 (x) + x]dx 

+ J S _ Oie^s^dsdz 

™ (A £ - A ) j s _ ldy + £ p e (8)d9f:l L {V 2 (x) + x]dx + 0(\ ln £ |e§), 

(3.98) 

as e ->• 0. Similarly, recalling JOS]), (|3~T5| . (IX57)) . and (|3~55)) . we have 

(u* - A(eiyj) dy = (A E -A ) / ldy + / " /3 £ (0)d0 / [V 2 (x) + x]dx+0{ei), 

V ' JS JO J-2L 

(3.99) 

and 

/ u 2 £ dy = [ * /3 £ (8)d8 [ V 2 (x)dx + 0(ei), (3.100) 
is + Jo Jo 

as e — )• 0. Moreover, thanks to (|3.90[) . (|3.91l) . we have 

(u 2 .-A( S %y))dy = (\ £ -\ )[ ldy + 0(e»), (3.101) 



■D e \(S US_) x ' JV E \(S US-) 



and 



lR 2 \(V e uS + ) 

as e — > 0. Let us keep in mind that 



/ u 2 e dy = 0{e- ce § ), (3.102) 



ui{y)dy = e-* / ^(y)dy = e"* . (3.103) 

JR 2 

Furthermore, recalling (|1.6I) . (|1.7[) . (|1.9[) . and (|1.25[) . for small e, we can write 

/ A+(y)dy= [ A+(y)dy- [ A+(y)dy+[ A+(y)dy, 
Jv E Jv Ju e Jv e 

where U £ C P , V £ n X> = 0, |W e | + |V e | < C\\ne\h, and A+(y) < C|lne|5 e if 
yeU £ UV £ . Thus, via (fL5]> . (pH]) . we infer that 



/ A + (eSy)dy = E - ^+0(|lne|ef) ase^O. (3.104) 
Jix 
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By combining (|3l98]) - (|3.104p . we deduce that 

(X e -X )\V s \+^ e /3 £ (e)d9j (^J [V 2 (x)+x]dx + V 2 {x)dx^j = C(|Ine|e») 

as e — > 0. Now, the validity of estimate (|1.45[) follows readily by noting that the 
sum of the above two integral involving V is zero. This can be seen by multiplying 
(|3.15p by V x , integrating the resulting identity by parts over (— oo, 0) and (0,oo) 
respectively, and recalling p,18|) . 

To finish, utilizing all the above, we will establish the validity of estimate (|1.46p 
for the energy of r] e . It is straightforward to see that 

Ge(Ve) = J v |i|V W | 2 + + L±W{eiy)u^dy, 

where u(y) = rj s {eiy) is the solutions of the stretched problem (|3.ip . Motivated 
from (jl.lOp . and recalling (11.51) . it is easy to check that we can rewrite the above 
relation as 

Ge(Ve) = U M2 \Vu\ 2 dy+^fv e (u 2 -a(eh))\y-±(J M M + ) 2 dy)e- 2 + ^e- 2 

(3.105) 

Similarly to the above proof of (ll.45p . keeping in mind the proof of Proposition l3.61 
(H~45l) . ([3~18l) . (I3T47)) . ([3~49]) . (ET87)) . and relations gl}, gl]) below (whose proofs 
do not require (|1.46|) ). we get 

J R2 \Vu\ 2 dy = f s _ (et/3 4 T/ 2 (/3.s) + 0(e*)) (1 + eiks)dsdz + 0(1) 

= s:™i v 2 ( X )dx+o(i) 

as e — > 0. Furthermore, as in the above proof of (|1.45p . we have 
^3 J_ (y?- a (eh)) 2 dy = 0(l), 

and 

-- r 2 — - r 

e —r- / (u 2 -a{eiy)) dy- £ —^ a 2 ( £ iy)dy = 0( e S), 

as e -> 0. Now, the validity of (|1.46p follows at once from (|1.45l) . (I3.105p . and the 
above three relations. Alternatively, we could have used the formula 

G e (r, s ) = ^ s -±J^dy, 

which follows easily by testing equation (|1.8p with r] e . 
The proof of Theorem II .11 is complete. □ 
We now outline a few remarks. 
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Remark 3.9. The maximum principle yields the upper bound: 

Ve(y) < max y G M 2 , 

(see also iS.lOb}) below). 

Remark 3.10. As in 125 ], where the authors refer to an idea of Shafrir, we can 
rewrite U.8\) in the form 

-e 2 A( y /a £ ~- rje) + n £ (q £ + ^/a^){r] E - ^fa~ £ ) = -e 2 A( v /a7) in V e . 

The above relation suggests the following, which can be proven similar to |150] : We 
have 



»fe(y) = v^ + g o +o(0. ( 3 - 106 ) 

where e _2 o(e 2 ) — > uniformly on any compact subset ofT>o as e — > 0. Keeping in 
mind Ji.7| ) which implies that 

„A(\/a7) / o, , 5\ 
£ 2 w = ( £ 2 f - 2 uniformly in Vq as e -> 0, 
2a £ V / 

and Ji.^flp , we are tempted to believe that hS.lOb}) can be extended to hold uniformly 
in the domain T> £ \{—Kes < t < 0}, with K large, if e — > 0. A possible approach 
could be by seeking a more refined inner solution with V + s4> in place ofV in i3.30\) . 



where <p is determined by solving a linear equation of the form M.{4>) — f(x, z), x € 
K, z € [0,£~3^ e ) with f known (in terms of the curvature k £ , V, a £ and their 
derivatives) and A4 as in H3.20)) : we refer to [135 for a related problem. Actually, 
we have computed that in the radially symmetric case, in N > 1 dimensions, we 
have 

M{<j>) = ^r 1 ^ - \a„{R e )[3- A x 2 V, 

where R £ is the radius of T> £ and (3 — [—a r (R £ )]i. Moreover, due to matching 
conditions with a £ {R £ + e$ ft^x), we need that 

1 3 

<j>{x) + -a rr (R £ )(3^ i (~x)^ ->• 0, x -> — oo; 0^0, x -> oo. 

In the special case of the model harmonic potential, this lower order term in the 
inner solution has been formally derived in [94] and rigorously in 104 . It might 
also be useful for the reader to take a look at j4-%4\ ) below. 

Remark 3.11. As in |125[ Prop. 2.1 e)] , it follows that 

\\r) e — v^ellc 1 ^) ^ C^e 2 f or an V compact subset K. C T>q, 

if e is small. 

Remark 3.12. Given A > mingN W , problem (|1.22[1 . in N = 2 dimensions, has 

a unique positive solution, for small e, as has been proven recently in [125] by com- 
bining ideas of Brezis and Oswald [43j ( see also the uniqueness part of Proposition 
\B.S\ herein) with those used in the proof of De Giorgi's conjecture in low dimensions 
[107] . This fact allows us to work exclusively with equation \1.8\) . since positive 
solutions of the latter coincide with the unique real valued minimizer of G £ in H. 
On the other hand, using \1.27\ ), which holds for every positive solution of (|1.22[) _ . 
it is easy to see that the method of [125J can be extended to the case of arbitrary 
N > 1 dimensions. Hence, our Theorem \1.1\ can be extended naturally to treat the 
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case where the functional G £ is considered in arbitrary dimensions, with the analo- 
gous conditions on the potential. In the radially symmetric case, uniqueness results 
which allow the case where X = inf K w W may be found in [3] and |104j . 

Remark 3.13. One can also prove an analogous result to Theorem \l.l\ for the real 
valued minimizer of G e in 

3 = |uG <' 2 (X> ;C) : jf M 2 dy=lj. 

This problem has been studied in [130], in a three-dimensional setting, with potentials 
of the form hi. 10) . The special property that AW > in T>q was used in an essential 
way in the latter reference for estimating the minimizer near the surface dT>o, along 
which it has a steep corner layer. We also refer the interested reader to [8] and [13] 
for the case where T>q has annular shape. For a numerical treatment of the problem 
we refer to [29] . As in Subsection \3.1\ below, it is not hard to see that in this case 
the layer profile near dT> should be determined by the unique solution of 

v X x - v{v 2 + x) — 0, X < 0, 

(3.107) 

v{x) — \J — x — > as x — > — oo; v(0) = 0. 

We refer to AvvendixW\ below for a treatment of the above problem in relation with 

The minimization of the functional G e in W ' (Vq), subject to the mass con- 
straint, leads to the equation in (|1.22|) _ with Neumann boundary conditions. The 
latter singular perturbation problem may be treated by using in place ofV, in i3.30\) . 
the (reflection of the) solution described in Remark \B.6\ below. 

Remark 3.14. Identical estimates to tl.38\) - [l.42\ ) (with t replaced by the signed 

distance from dT>o) hold for the solution fj E of (j 1 . 22[1 (with X — Xo, q — i). One 

can use the function J 1 ' € % as a competitor in order to give a self-contained 

WVe UL 2 (K 2 ) 

proof of fZHiP ) from which il.2.5\) follows readily (see jlOj . [125j ). (As in the proof 
of |J-46p , keeping in mind that ||%|| 2 : 2( R 2) — > ||^4 + ||l 1 (r 2 ) 7 the main contribution 
would be from the gradient term). 

Remark 3.15. In the case where the potential W is of harmonic type, as in (1. 16]) , 
it was observed by the authors of [125] that 

i v VAo + X e _ / v%y \ _ A e 

Ve{y) = 7= Vi 7t ; , with e = 



Ao VvAo + X e J Ao + X e 

where f)g was defined in Remark \3. 14\ This identity and a technique of Struwe [1251 
Lemma 2.3] were used essentially in their proof of |1.^5[ j for this special class of 
potentials. 

Remark 3.16. Relation ljl.44\ l implies that, if e is small, the linearized operator 
C e about u e is invertible. Hence, the implicit function theorem (see for instance 
[19j ) implies that there exists a small number Eq > such that, in addition to being 
isolated (for each e), the minimizers 77(e) = rj e depend smoothly on e € (Q,Eq) (in 
all the usual function spaces). In particular, we have 77 : (0, £0) — > W 1,2 (M. 2 ) is C . 
This last property yields at once the first part of Lemma 2.3 in 125 , mentioned in 
Remark \3.15\ above. 
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Remark 3.17. In the special case of the model harmonic potential W(y) — \y\ 2 , 
since r\ e is radially symmetric [108 ], we can define v e : (— oo,A e £~3] — > M by 

Letting £ = (A e — |y| 2 ) /e^ , then the equation in U.8\) becomes equivalent to 

4(A e - ei^d^e ~ 4e%d 6 v e +& e - v\ = 0, f e (-oo, A e e _ ^]. 
Af /irsi glance this might look rather counterintuitive but this strategy, already used 
in [104] . allows one to use directly (2A e )sy ^— - — as a global approximation 
(recall fOSU - fQ3l) ). 

Remark 3.18. In £/ie radially symmetric case, estimate l{1.4(ty was formally pre- 
dicted in [67, 65] , and rigorously proven and extended very recently, for the case of 
the harmonic potential, in [105] at the same time that the current paper was written. 

4. Further properties of the ground state r) e 

As a byproduct of our construction of the ground state r\ e , we can extend to the 
non-radial case relation (|1.30[) . improve relation (| 1 . 29[) . show that r\ e has maximal 
Holder regularity, and improve relation (jl.28l) . Finally, under an additional but 
natural non-degeneracy assumption on the potential, we can refine bound (|1.39[) . 

Corrolarry 4.1. There exist small constants c, d' > such that, given D > 0, we 
have 

(rie)t<-c(\t\+e%)-i (4.1) 
in i—d! <t < Dei , 9 € [0, i E )\, provided e is sufficiently small. 

Proof. We will present it for the solution u e of the stretched problem ([3.1)1 . It 
follows from ([3.79)1 that, if e is small, the function tp* — u £ — u ap satisfies 

r Ce(\s\ + if - 10fc-i <(3s<0, 

\Ay<P*\<{ (4-2) 
[ Cee~ c ^ if < [3s < 10fc"i. 

In establishing the above estimate, the term that needed extra care was 

which can be estimated by noting that 

u 2 ap +\a{eiy)\ < Ce§(|s|+1) and \<p*\ < Ce(\s\ if {\f3s\ < lOSeS}, (4.3) 

(recall (IPS)) . (|337[) . (|335[) . and ([37T7[)l Making (mild) use of gTJ and the second 
estimate in (j4.3[) . via standard interior elliptic regularity estimates )109j (applied on 
balls of radius one) or the interpolation- type inequality of Lemma A.l in [38], the 
statement of which is included below as Lemma 14.11 for the reader's convenience, 
we obtain that 

\V v <p*\ < Ce (4.4) 

in the neighborhood of f e described by {\j3s\ < 5Se~ i}, provided e is sufficiently 
small. 
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By (|3T7l) . (pHU)) , (j3~49l) . ([3~51~j) . and gl}, we deduce that, given D > 0, 

(u e ) s = e^/3 2 V;(/3 S ) +0(e) < -ce*, (4.5) 

uniformly in the neighborhood of r e described by {— L < x < D}, as e — > 0. In the 
same manner, recalling (|3.19l) and the second estimate in (|3.40[) . we find that 

(u e ) s = e^/3 2 K(/3 S )+0(e| S |5) 

< -ce*\s\-^ +Ce\s\^ (4.6) 

< -|£3|s|-5, 

as long as — de~i < (3s < —L, for some small constant J, provided e is sufficiently 
small. The corresponding assertion of the corollary, for the solution of the equivalent 
stretched problem p. II) . follows readily from (I4.5[) and (|4.6[) . 

The proof of the corollary is complete. □ 

The following is Lemma A.l in 38 : 

Lemma 4.1. Assume that u satisfies 

-Au = f in SJcl". 

Then 

|v u (2/)| 2 < |||/IU~<n)IMU~(n) + * anJ Hi-C")} v ye n - 

where C is some constant depending only on N . 

We can also improve the bound (fl~29|) : 
Corrolarry 4.2. If e > is sufficiently small, we have 

\\Vr)e\\L°°(M*) < Ce~%. 

Proof. We will prove the corresponding assertion for the solution u e of the stretched 
problem ([34]). From the proof of Proposition ELS recalling ([3~19| . (|347|) . ([348)) . 
(j3"4T)l) . and (|4"4) . we find that 

|Vu(y)|<C E 5 ifye{|a:| <3fe-i, ze [0,e-*4)}, (4.7) 
if e is small. Furthermore, recalling (|3.77[) . we get 

„.2 



a(eiy) =2(p*yja(e*y) + <pl = 0(e 2 ), 



uniformly in T> e \{— 25e 3 < x < 0}, as e — > 0. Hence, by the equation (|3.1[) . we 
obtain that 

|Ati| < Ce^ in £> e \{-2fe-i < x < 0}. 
Consequently, by the interpolation- type inequality of Lemma A.l in |38j (see Lemma 
14.11 above ), we infer that 

|Vtt| < Cei in V e \{-36s~% < x < 0}. (4.8) 
Similarly, recalling (|1.4j) and ()1.41|) . we have 

|Vit| < Ce~ ce ^ outside of V e U {0 < x < SSe'i}. (4.9) 
The corresponding assertion of the corollary for u follows readily from (|4.7|) - (|4.9j) . 



The proof of the corollary is complete. □ 
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In the following corollary, we will show that r\ e has the maximal Holder regularity 
available (recall Definition Q] from Subsection II. 2p . 

Corrolarry 4.3. If e is sufficiently small, we have 

||»kllov>(R») < a (4.10) 

Proof. From (OS]) . (|l~39|) . (|3~18)) . (j3l9l) . (j3~49)) . and abusing notation, it 

follows that 

' c(\t\+eiy < % <c(\t\+eiy , 

\(Ve)t\<c(\t\+e^y\ (4.H) 

\(Ve)e\<c(\t\+ei)\ 

in the region described by {\t\ < 35}, if e is small, having further decreased the 
value of (5 if necessary. Abusing notation once more, let y, = (ti,9i), i = 1,2, with 
\ti\ < 35 and Qi <G [0,£ e ), be any two points in that region. We write 

- %(y 2 ) = A) - Ve(t2,0l) + Ve{t2,0l) ~ Ve(t2,e 2 ). (4.12) 

Now, instead of considering the difference % (i i , #i ) — ?y e (i 2 , $i ) , we will first consider 
the difference rjl(ti,9i) — 77^ (£ 2 , ^i)- We have 

r£(*i,0i)-»£(t2,0i) = (ti-*a) / 2 Wt (i 1 +r(< 2 -i 1 ),# 1 )dr. 

Jo 

So, thanks to (l4~TTj) . we find that 

|»£(ti,0i)-»£(ta,0i)| <C|ti-ta|. 
In turn, via the lower bound in (|4. 1 1|) . the above relation yields that 

\r,Si,6x) - VefaM < C ^r^'j ti - *al* < - <a| + |0i - 0a|)* • 

+ I = 

(4.13) 

Similarly, we obtain that 

|^(t 2 ,0i)-7? e 2 (*2,02)| <c(|* 2 |+e§) Iflx-flal, 
which, as before, implies that 

1 

\Ve(t 2 ,0l) ~ Ve(t 2 ,6 2 )\ < C (\t 2 \ + £§) 5 |6» x - 2 1 < C (|*i - fc| + |0i - 2 |)' ■ 

(4.14) 

Hence, by (|4.12[) . (|4.13p . and (|4.14l) (also keeping in mind Remark (|3 ^) . we deduce 
that 

if e is small. That was the hard part. In the remaining regions of the plane, by 
virtue of (|4.8j) and (|4.9j) . we see that |Vry e | < C which implies that, in those regions, 
the family r\ e is in fact uniformly Lipschitz continuous, as e — > 0. 

The proof of the corollary is complete. □ 
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Remark 4.1. Uniform Holder C°' Q , < a < 1, bounds for Gross- Pitaevskii sys- 
tems, where the singular limit functions have Lipschitz regularity (a = 1 ), have been 
proven in }1671 by blow-up techniques and the monotonicity formulae of Almgren 
and Alt, Caffarelli, and Friedman (see also [53j ). The importance of our result lies 
in the fact that the Holder exponent 1/2 in (14.101) equals to the exact maximal 
Holder regularity of the singular limit profile. To the best of our knowledge, this 
property has been proven, in singular perturbation problems, only in one- dimensional 
problems, see . 

The following corollary answers a question posed to one of us by A. Tertikas in 
relation with |191j . 

Corrolarry 4.4. Given a G [0, |), we have 

r] E ^y/A+ in C a (R 2 ) as e ->■ 0, (4.15) 
but r\ e does not converge to V ' A + in C^(R 2 ) as e — > 0. 

Proof. Given a £ [0, h), in view of (|1.13p . Corollary |4.3l and the compactness of the 
embedding C^(2V ) ^ C a (2V ) (see ^Q9\), we find that r? e -> y/A+ in C a (2T> ) 



as e — > 0. Now, the desired relation (|4.15l) follows via (|1.42l) and 

On the other hand, the following simple argument shows that we do not have 
convergence in C?. Let {s n } be a decreasing sequence such that e n — > 0, and, 

2 

abusing notation, consider the points y n = (t n ,9 n ) = (— £n,0). If n < m, thanks 
to (TOS) . (fT35jl . and JO}, we get 

1 

I, _ || he„(y n ) - ^,„(y m )| At - ej V § 

I'd \c n ~r £m / 

where c, C are independent of n, m. Now, choosing for example e n — (^) 3 and 
m = 2n, we conclude that is not Cauchy in C'. 

The proof of the corollary is complete. □ 

Remark 4.2. VFe expect that there exists some constant C > swc/i that, for small 
e > 0, we have 

\\ Ve -VA+\\c^)<Ce^- a \ 0<a<l. 



2' 

see also Remark \4.4\ below. 

The following result is motivated from (|1.28l) : 

Corrolarry 4.5. Given < a < 1, if £ is sufficiently small, and A as in U.9\) 

there exists a constant C > such that 



\t]e - V A+| < C + e _ * Q |A e - AolJ < C*(e^ Q + £ 3-3«)Vyl+, (4.16) 
ai points in Dq whose distance from dT>o is greater than £§ a . 

Proo/. From ([O]). ([T7j) . (|T2l)]) . and P^lDj) . it follows that in the region 
£> \{-^ Q < * < 0} we have: 



V^y)>ce*, (4.17) 



y/Wi- - W(y) = 0(\X e - A |e-«), (4.18) 
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and 

Ve(y) - VA £ - W(y) = 0{e 2 -i a ), (4.19) 
uniformly, as e — > 0. Now, the assertion of the corollary follows readily by combining 
(|1.45j) . Remark ll.31 and the above three relations. 

The proof is complete. □ 

Remark 4.3. Note that estimate fi4.16\ ), when a — \, considerably improves esti- 
mate U.28\) . which was originally proven in [8] (for solutions of il.8\) with X £ = \q ) 
and in the sequel used in [10], [13], and |125) . 

Remark 4.4. By hl.38}) - fl.39\) , and 14-16] ) with a = 1, for small e > 0, we obtain 
that 

\\v e - VA+\\ L ~pi) < Cei. 

Under an additional but natural non-degeneracy condition, satisfied by most 
potentials used in physical applications (recall the discussion in Subsection ll.2[) . we 
can improve estimate (|1.39[) . 

Proposition 4.1. If we assume that W € C 2 and 

W 4t (O,0)>c>O, 6e[Q,£), (4.20) 

then 



r le {y)=E*(3 e {6)V ( 



i+o( E t)(4 



(4.21) 

uniformly in {0 < t < d, 8 £ [0,^ e )}, as e — !• 0, where d > is some small constant. 

Proof. Once again, we will work with the equivalent problem in stretched variables. 
Our aim is to estimate ip* = u e — u ap using equation (|3.79|) , as we did for estimate 
(j3~78| . By virtue of ([3T5]) . (j3^3| . (|3T24]) . (p33]) . and (j33T|) . for small e, we get 
that the remainder in (|3.80|) satisfies 

\E\ < Cex 2 M if x G [l,2fc~i). (4.22) 

In view of (jl.32l) . (|3.5T[) . and (|4.20p . it is easy to see that, decreasing 5 if necessary, 
we have 

(Zulp - a (£h)) > P 2 x + ceix 2 if x G [1, 25e~i), (4.23) 

provided e is sufficiently small. We point out that in the above relation the constant 
c is independent of small 5, s. Let > be determined from 

- f" + x^ = x 2 M- *(1) = 1, 4>(oo) = 0. (4.24) 

In order to proceed, we need some estimates for ^ . A short calculation shows that 

, s , rt 2 {A\fdt 

* = (Ai /i with ti = Jx /t \, — . 

(Ai) 2 

Note that, from (|3~22|) . (f3T24|) . we have 

Ai' ~ -x^ Ai as x -> oo. (4-25) 

By the way, a neat way to show the above relation is to use L'hospital's rule to find 
that 

x-HAQ 2 
hm — — — - — = 1. 
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Now, refereing to L'hospital's rule once more, we get 

h' ~ —x? and h ~ —x? as x — > oo. 
2 5 

Hence, we find that 

* — -x%Ai and ~ -ia; 3 Ai as x -> oo. (4.26) 
5 5 

Keeping in mind that x = /3(etz)s, via formulas (|3.6|) - ([3.7|) . we get that 

A y ty = /3 2 *" + <D{e^)x^>' + C(e^)x 2 1'" + 0(e*)V', (4.27) 

uniformly in jl < x < 25s~i , z € [0,^ e e _ ^)|, as e — > (note that here O(-) is 

bounded uniformly in small 5). Using ([177). ([3760]) . (|4~23|) - ((47271) . and further 
decreasing 5, we readily find that 

-£ e (*) > /3 2 x 2 Ai + c£ix 2 * + 0(et) x 3^ + ( e |) x 4 Ai + C)( - £ §> )a ,3 Ai 

(4.28) 

> ^-x 2 Ai + ceix 2 * 

if a; € [l,2fc~i), provided £ is sufficiently small (c, £>(•) independent of small 5). 
It is nice to note that, in the above calculation, the cubic power in the second 
asymptotic relation of (I4.26[) was "the most appropriate" one in order to absorb 
the last term of (|4.27|) into the term /3 2 x 2 Ai, by decreasing S. Similarly, keeping in 
mind (|3722]) - ([3723]) , we find that the function 

B(y) = Bi(x-2fe - i), ye {x E [l,28e~*), ze [0,e~h e )}, (4.29) 

satisfies 

-C E {B) > -/3 2 B" + f3 2 xB + ceix 2 B + O0)x 2 {x-26e-i)B + O(ei)B' 

> 25(i 2 e-^B + ce^x 2 B + 5 3 0(e-^)B + 0(e*) (\x\v + Sh~^ B 

> Sp 2 e~iB + ce%x 2 B, 

(4.30) 

if e is sufficiently small, having further decreased S if necessary. From now on we 
will fix 5. In view of (|37%T|) . (0722]), (0725]), and (gTSTJ, given M > 1, the function 

$(y) = Me |*(x) + Bi(x - 2fe"*)| 

satisfies 

-£ £ ($) + 7V($) + S > M^£a; 2 Ai + cM£la; 2 4' + M ( 5/3 2 £3B + cM£fa; 2 S 

-C*£a; 2 Ai - CM 2 £§x 5 t/(Ai) 2 - CM 2 e*x*VMB 
-CNPe^VB 2 - CM 3 £3a;iAi - CAPe^B 

> M^ex 2 Ai + M5p 2 e^B + aMei x^Ai ~ C 2 M 2 e^x 5 V(A 

for some constants C\,C?. > (independent of e), if £ < e(M) is sufficiently small. 
Note that the seventh, eighth, and tenth term in the first inequality's righthand 
side were absorbed into the corresponding third term, whereas the ninth into the 
first. What we want to do next is to somehow "get rid" of the last term in the 
above relation, and end up with a positive righthand side. We will achieve this by 
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absorbing that term into the one that proceeds it. By virtue of (|3 . 23[> . p.24[) . and 
(|4.25p . if M is sufficiently large, there exists xm > such that 

C 2 Mx?V Ai < ci if x > x M - (4.31) 
(Note that xm — > oo as M — > oo). Hence, it follows that 

- £ e ($) + AT($) + £ > M5f3 2 e?B + M-^ex 2 Ai > (4.32) 
in the strip-like domain described by 

S £ = [x M < x < 2fe"3, ze[0,4£^)}. 
Now, in view of (j3~75j) . (g^BJ), and (|43T|) . we can fix a large M > such that 

ip* < $ on {a; = xm}; ^* < <& on {x = 2fe~^}, (4.33) 

if e is small. By (|4.32|) . ()4.33|1 . and a standard maximum principle argument, making 
use of (|3.56p and the property that 

\N($)-N(<p m )\<CE*\$-tp m \ mS £7 

we deduce that tp* < $ in S ei if e is small. Similarly we can show that tp+ < — $ 
in <S e , if e is small. The desired assertion of the proposition (for the equivalent 
stretched problem) now follows via (jl.39l) . Q3.23p . (|4.26p . and noting that 

2 EH 

Bi(x - 25e~*) < Bi(-5e~*) < 2Ai(fe") < 2Ai(x) 

if x G [xmi 8e^i) with e sufficiently small. 

The proof of the proposition is complete. □ 

Remark 4.5. Note that the first term in the righthand side of dominates 
for < t < el. 

Assuming additionally that W is radial and convex (outside of T> e ), we can derive 
an explicit global upper bound on the minimizer. 

Proposition 4.2. Assume that the potential trap W is radially symmetric with 

W rr {R E ) > 0, and W rr (r) > if r > R e , 
where R £ denotes the radius ofT> £ . Then, we have 

Ai ( fi £ s ~^' \ 

Ve(s) < / £ e5 ( ve(r) V s>r>R e . (4.34) 



In particular, it holds that 

»fc(*)<e*(& + o(l)M*(&^^J ifs-R e ^ei, (4.35) 

as e —> 0. 
Proof. Since 

-a e {r) = W(r) - W(R e ) > W r (R £ ){r - R £ ), r > R s , 
we see that the minimizer rj e is a positive lower-solution of the linear equation 

" e% [Vrr + Ivr) + - R e ) V = 0, (4.36) 
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if r > R e (recall that (3 S = [W r (R e ))^ > 0). On the other side, making use of 

(|3.22p and the fact that (Ai)' < 0, we readily find that Ai (Pe 1 —?^ is a positive 

upper-solution of (|4.36|) if r > R £ . Hence, by the maximum principle, we deduce 
that relation (|4T34|) holds true. In turn, via pT39]) . ((3~23l . (|3~2"6]) . and Proposition 
14.11 relation (|4.34p implies the validity of relation (|4.35[) . 

The proof of the proposition is complete. □ 

Remark 4.6. It seems plausible that the techniques of the very recent paper [59] 
can be extended to derive a WKB (Wentzel- Kramers- Brillouin) type estimate, in the 
region (R,oo) (with the obvious notation), for the ground state solution of U.2S\) 
with q = 2 and N = 1 . 

5. Refined estimates for the auxiliary functions £ £) f £ in the case of 

RADIAL SYMMETRY 

In this section, restricting ourselves to radial potentials with T>o a ball, building 
on our previous results for the ground state r\ E , we will improve upon the estimates 
obtained recently in [10] for the auxiliary function f e in (|1.31[) . As we have already 
discussed in Subsection 11.31 the latter estimates were essential for the analysis of 
[TO] regarding the functional E e , defined in (jl.151) . We believe that the improved 
estimates herein may provide important intuition for the treatment of the general 
case, which may ultimately lead to the resolution of the open problem raised in [10] 
(recall the discussion in Subsection II. 3j) . 

In the general case, for potentials as described in Subsection ll.il we define £ e to 
be the solution of 

div (4 VC ) = ~ 2 ' y G ^ (y) ~* °' lyl ~* M ' (5 ' 1} 
so that V £ e = x T/g. An integration by parts in (]1.19|) yields 



FeH = 141 U Vw|2 "T JW J + & (H2 " 1)2 i dy ' (5 - 2) 

where Jw — x (iw, Vw) = (iw yi ,w y2 ) is the Jacobian. 

We recall that the function f e :— t; £ /Ve' appearing in the functional F s , is im- 
portant since it is well known that vortices in the interior of T>o first appear near 
where this function attains a local maximum 8, 9, 1251 1126] : its importance is also 
clear from (]5.2|) . since it controls the relative strength of the positive and negative 
contributions to F e . 

In the case where the potential W is radially symmetric, one can solve problem 
([5.1]) explicitly to find that the functions £ e , f £ are given by relation (jl.3ip . In 
particular, if the domain T>q is a ball, it has been shown in [10] that, for small e, 
we have 

Cdist(y, dV a ) + Ce% if y € V , 
fj\y\) { and ||/ E -/ |U-(R) < Cei, (5.3) 

Cea if not, 



where fo is the function in (|5.7p below, which solves the "limiting" problem corre- 
sponding to ([5.ip : 

div - -2 in V , - on dV , (5.4) 
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with A as in (|1.9[) . Existence and properties of a positive solution £q of (15.41) have 
been established in [12] , 

The following proposition refines and improves relation (|5.3[) . 

Proposition 5.1. // the potential W is radially symmetric with Dq = {y € K 2 : 
r = \y\ < R}, then the function f e , defined in fl.31)) , satisfies 

f e (r) = R^-h^V- 2 (p/—^) r V 2 (a)da + o( e i), (5.5) 

uniformly in [R — o(ei), oo), as e — > 0, where R £ denotes the radius of the ball D e . 
satisfying H5.8\) below, and f3 e = [W'(R e )]3 (recall U.36\) ). 
Moreover, if e is small, it holds that 

||/e-/o|U-(R)<CE*, (5.6) 

where 

( if r R sA( S )ds, 0<r<R, 
fo(r) = I (5.7) 
{ 0, r>R. 
(In view of i U.6'1) . and \1. ty , an application of L 'hospital's rule shows that /q(_R _ ) = 

2 /' 

Proof. First of all note that, thanks to (|1 .45[) . we have 

R E = R + Q{ei) ase-^0. (5.8) 

By virtue of (|1.4[) . (|1.6p . (11.71) . and (|1.45l) . there exists a constant c > such 
that, given K ^> 1, we have 

r)l + W{r) - X e > W(r) - \ E > c{r - R) 2 +cKe%, r > R + Ks%, 

provided that e is sufficiently small. (We note that K is considered fixed in the 
corresponding relation in |10j). Then, by a standard barrier argument in equation 
(|1.8[) . we deduce that 

Ve(s) <Tie{r)e- KiE ^ {s2 ' r2) , s>r>R + Ksi, 
if K is sufficiently large and e sufficiently small. As a result, we get 

f £ (r)< se- 2Kle ~^^ 2 - r ^ds=jK-^e^,r>R + Kei, (5.9) 

J r 

if e is small. 

If r € [R — d', R + Kez], in view of Corollary 14.11 and (15.91) . for small e, we have 
Mr) = ^^hn 2 Mds+ r ^^f e {R + Ke^) 

(5.10) 

= ^J r R+Kj sn 2 (s)ds + 0(K-ish 

uniformly as e — > 0. If r e [R, R + Kss], from ([09]) . (|3T8]l , and (JSHJ), it follows 
readily that 

e3 (5.11) 

+0 ( £ §y- 2 (^) /~ y 2 (a)da) + Ojf (eS), 



THE GROUND STATE OF A GROSS-PITAEVSKII ENERGY 



47 



uniformly, as e — > (the constant C_fs"(l) may diverge as K — > oo). The second term 
in the righthand side of the above relation can be estimated as before, by noting 
that, thanks to (]3.15|1 . we find that 

V(a) < V(K)e- Ki{,J ~ K} 1 a > K if K is large. 

Thus, by l|5.10l) . (15 . 1 1[) . we infer that relation (|5.5[) holds true in [R,R + Ke?\. In 
fact, by (|5.9j) and the above relation, we deduce that (|5.5[) holds true in [R+Kez , oo) 
as well. If r € [R — d, R], similarly as before, but this time using (|1.38[) instead of 
(11.391) . we arrive at 

f £ (r) = R e ^eiV-*(p^)S™^V*(*)da 

(5.12) 

+0 (\r - R £ \ 2 + e\r - R £ \i + ei\r - R e \^) +o(ei), 



uniformly, as e — > 0. The above relation implies at once the validity of (|5.5j) in 
[R, — o(e3),R] as e — > 0. Consequently, we have established the validity of f|5 . 5|) . 

Next, we will show the validity of estimate (|5 .€>[) . If r 6 [R — e~s~ , R], with 
I < a < 1, recalling (fl~7)) . we have c(i? - r) < A(r) < C(R - r). So, as in [10], we 
obtain that 

C [ R 

Mr)< / s(R- s)ds <C(R~r) <Ce^. (5.13) 

R - r J r 

Furthermore, thanks to (|3.18[) , (|5.5p . if e is small, we find that 

f s {r)<Ce^, r G [R-e^,R]. (5.14) 
If r € [0, R- £-# ], following [10], we write 



(5.15) 



/ B (r)-/o(r) = sv^ds-^Sr sA( s )ds 

= I + 11- III. 

Using (|1.6[) , (|1.7|> . (|1.40p . Corollary 14.11 and our earlier estimates on f E , f for 
r > R — £ 3 j we see that 

|//| <Cf e {R-e*r) <Ce^ and |JJI| < C/ (i? - e*) < Ce*. 
We further decompose the remaining term as 

'-(wrm) C + W) ST" ' Ws) - m d '~ 

Using Corollary 14.51 if | < a < 1, for small e, it follows that 

|/| < Ce 2 ~ 2a ( R ~ e 3 s^-ds + Ce 2 - 2a ' s^fi-ds. 

Due to CorollaryE] we have 2fl4 < 1 if R-d' <r <s< R-e^ . If < r < R-d', 

2/ \ E 

then rjl(r) > c, and so ^f-^j < C. Thus, the first integral in the above relation is 
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bounded by Ce 2 2a . The second integral is estimated similarly, using (11.71) instead 
of Corollary 14. II Therefore, relation (|5.15|) implies that 

\fe(r) - fo(r)\ < Ce 2 ~ 2a , r G [0,R- 

provided that e is sufficiently small. The validity of estimate (|5.6I) follows at once 
by combining (15 . 13[) . (|5.14l) . the above relation, and choosing a = |. 

The proof of the proposition is complete. □ 

6. Open problems and future directions 

What follows is a list of questions which are currently unresolved. These are 
presented as an illustration of where our interests lie. No attempt is being made to 
be precise in their formulation. 

A question that comes naturally to mind is to examine whether the estimates of 
Theorem ll.il for the minimizer rj e of G e inH, can be used to answer the interesting 
open problem posed recently in [10] . As we have already mentioned in Subsection 
11.31 the latter is to see to what extend the analysis of [TO], for the functional E £ 
in (|1.15p . continues to hold if one drops the assumption of radial symmetry on 
the potential W. Hopefully, our estimates for r\ e can be used in estimating the 
corresponding auxiliary functions £ E , / £ , arising in the functional E £ as in (|5.2j) . 
which for the radial case were given by (jl.31[) . It seems that the elliptic problem 
(|5.1[) . which determines £ e , seems to be a singular perturbation problem of its own 
independent interest. 

In the special case of the model harmonic potential, an approximate solution for 
(|1.22[) _. "close" to y/(X — W) + , of arbitrary order accuracy was constructed in |104j 
(keep in mind Remark I3.17|) . We feel that it would be very interesting if one can 
do the same thing for the case of general potential. A major difficulty (or problem) 
is that each term of the inner expansion diverges polynomially in a complicated 
manner as the distance from T increases (recall Remark l3.10p . see also the appendix 
in 54 . The construction of arbitrary order approximations is especially important 
in the treatment of singularly perturbed elliptic problems involving resonance, where 
the order of accuracy of the approximation is dictated by the space dimension, see 
for instance [160) . Problems of these type which feature the presence of a corner 
layer (similar to the problem at hand) have been studied recently in |135j (in two 
dimensions), see also Remark [114] herein. It should be noticed that in Allen-Cahn 
or (focusing) Schrodinger type equations, where its possible to construct arbitrary 
order approximations (see [160 and the references therein), the phenomenon is 
exponentially localized, i.e, the corresponding terms approach certain constants 
exponentially fast. 

Relation (jl.44l) implies that the spectrum, in L 2 (M. 2 ), of the operator L e , defined 
in (|1.43p . is bounded above by — ce^ , for some constant c > 0, as e — > 0. We expect 
that, making further use of the estimates of Theorem ll.il one can rigorously "link" 
the spectrum of L 6 to that of the one-dimensional "limit" operator Ai , defined in 
(|3.20p . as e — > (see also relations (|6.ip - (|6.2p below), and thus provide a valid 
asymptotic approximation for the eigenvalues of L 6 . In particular, the difference 
between the first two eigenvalues, called the fundamental gap, is of importance since 
it determines the rate at which positive solutions of the nonlinear heat equation, 
corresponding to Q1.22p _. approach the first eigenspace of L e (see [55], |118j . and 
especially [44 J. In the case where W is the harmonic potential, a rigorous con- 
nection between the spectrum of L £ and that of M. (see the discussion following 
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relation (|A.16|) below), as e — > 0, has been made recently in |104j (see also [134] for 
a related radially symmetric problem). A possible approach for the general case, 
where the potential is as in the present paper, could be by mixing techniques found 
in the aforementioned references with those developed in [TTJ [58] for the study of 
the spectrum of multi-dimensional Allen-Cahn and related phase-field operators for 
generic interfaces. One could even carry out an analogous program for the spectrum 
of the linearization of the defocusing nonlinear Schrodinger equation (ll.21[) _ at the 
corresponding ground state, recall (ll.47[) . The latter problem is often referred to as 
the Bogolyubov-de Gennes problem in the context of Bose-Einstein condensates, 
see [103U137] for recent studies specializing on the model harmonic potential. In the 
latter references, for reducing the complexity of the problem, the authors linearized 
at rjo (recall (11.24p ) instead of the ground state %. In this case, the linear operator 
defined by the lefthand side of the first equation in (|1.47p . with 770 in place of rj £ , 
has also been studied in |101j . in relation with |191j . 

Excited states are solutions of ljl.22|) _ with zero set inside the domain T>\ = 
{y : W(y) < A}. In the Thomas-Fermi limit, e — > 0, the Bose-Einstein con- 
densate is a nearly compact cloud, which may contain localized dips of the atomic 
density. The nearly compact cloud is modeled by the ground state of the defocusing 
nonlinear Schrodinger equation (ll.21[) _. whereas the localized dips are modeled by 
the excited states. In the one-dimensional case, with W the harmonic potential, 
excited states of (|1.22j) _ which are approximated, as e — > 0, by a product of the 
ground state and m dark solitons (localized waves of the defocusing NLS equation 
with nonzero boundary conditions at respective infinities, which after a re-scaling 
solve the one-dimensional (|1.49|1 ) were constructed in |174) by a finite-dimensional 
Lyapunov-Schmidt reduction (for the latter see for instance the book [20]). Loosely 
speaking, these solutions have a corner layer at the points corresponding to dT>\, 
and m (clustering) transition layers in (— C\ lne|e, C\ \ne\e), as s — > (see also 
[61]). Studies in the case of radial symmetry have been conducted in j!19| . We 
believe that, at least in two space dimensions, analogous excited states can still 
be constructed without any symmetry assumptions on the potential, by employing 
the estimates of Theorem 11.11 (in particular (| 1 .44(1 ) and the infinite dimensional 
Lyapunov-Schmidt reduction of [79] (see also [201p . In this context, the dimension 
N = 2 plays an important role for the solvability of a Toda system, periodic orbits of 
which determine, up to principal order, the location of m closed curves in T>\ where 
the excited state changes sign. These curves should collapse, as e — > 0, to a closed 
curve in T>\ that may be determined by the arguments in [86], [151] . [166 (if N = 1, 
the interfaces collapse at critical points of W). We expect that, in the case at hand, 
the reduction procedure is more delicate than [79] because the corresponding lin- 
ear operator L e has small eigenvalues (see also [71] for a related finite-dimensional 
reduction). If N — 1 or W is radial, one could also try to construct "high energy" 
excited states of (|1.22p _. having an increasing number of layers of order 1/e, as 
e — > 0, in the spirit of [HE], [SS] and the references therein. We remark that the 
result of [55] relied on ODE techniques, but it is expectable that a similar result 
could be proven for higher-dimensional problems. On the other hand, in the one- 
dimensional case, solutions of (|1.22p ^ bifurcating from the trivial branch have been 
studied in [146] . |187j . and |196j (see also |139j and the references therein). Let us 
make a formal connection between these two different types of solutions (layered 
and small amplitude respectively). Consider the one-dimensional case with poten- 
tial W having a global minimum which is attained at a unique point, say at y = 0, 
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and satisfying 

W(y) = W(0)+c\y\ a + o(\y\ a ) as y -> 0, 
for some constants a,c > 0. Arguing as in [1341 Prop. 3.25], it is not hard to 
establish that, given m € N, the first m eigenvalues of the linear operator 

-e 2 d yy + {W{y)-\)I, (6.1) 

which corresponds to the linearization of (|1.22[) about the trivial solution, are of 
the form 

W(0) - A + /i !: e^ft" + o (e^J as e 0, i = 1, • • • ,m, (6.2) 

where {fii} are the eigenvalues of the "limit" operator 

-d yy + c\y\ a I, 

(these exist by |121l Thm. 10.7], and /ii — > oo as i — > oo). In passing, we not 
that formulas (|6.2[) improve the corresponding lower bounds found in Theorem 
1.4 of |101j . Hence, we see that the number of negative eigenvalues (counting 
multiplicities), namely the Morse index of the trivial solution of (|1.22l) ^. (see |134] 
for the precise definition), diverges as e — > (recall that A > W(0)). From a variant 
of Weyl's asymptotic formula, see for example [3H pg. 521], it turns out that one 
has 

/li ~ ci a + 2 as i — > oo, 
for some constant c > 0. We expect that, by refining the above argument, one 
can prove that the Morse index of the trivial solution is of order greater than or 
equal to 1/e, as e — > (keep in mind that Landau's symbol in (|6.2[) may depend 
on m S> 1). On the other side, it seems plausible that the operator in (|6.1j) does 
not have any negative eigenvalues if e is sufficiently large (by Poincare's inequality, 
this is certainly true when considered in a fixed interval with Dirichlet boundary 
conditions). Consequently, since the eigenvalues are smooth functions of e (by virtue 
of their simplicity [60] Th. 3.1, p. 482], see also |160| h we expect that there exists 
a sequence {e^} with e\ > £2 > • • • > £j — > as i — > 00 such that, for each e = £j, 
zero is an eigenvalue of the linearized operator described in (|6.ip . This suggests 
that the aforementioned local bifurcation of solutions of ([1.22|) _ . from the trivial 
branch, takes place at each e = We further expect that, using global bifurcation 
techniques |178j (see also [2]), one can show that these solution branches reach, as 
e —> 0, the layered solutions of (|1.22p _ that we discussed previously. (We point out 
that solutions belonging to the i-th branch have exactly i — 1 zeros) . We note that 
analogous eigenvalues of the form (|6.2j) . with the obvious modifications, also exist 
in the multi-dimensional case, and existence of many solutions for the nonlinear 
problem may follow by adapting Theorem 10.22 in [2Tj . Moreover, in the "flat" 
case (motivated from a definition in [47], see also |147| ). where the potential W 
attains its minimum value over a domain £Iq , we expect that the multi-dimensional 
operator, corresponding to (|6.1|) . has eigenvalues of the form 

W(0) - X + ^{^e 2 + o(s 2 ) ase^O, 

where {/ii(fio)} are the Dirichlet eigenvalues of —A in f2o ■ We believe that the 
eigenfunctions associated to the above eigenvalues and the corresponding eigen- 
functions of the Laplacian share the topology of their level sets, as in [90]. In any 
case, motivated from results in |llj . we believe that the existence and multiplicity 
of solutions to (|1.22p _. is strongly associated to the number of negative eigenvalues 
(counting multiplicities) of the corresponding linearized operator about the trivial 
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solution. In the radially symmetric case, the topological approaches of [111) . |187j . 
for constructing nodal standing wave solutions of the focusing NLS, should also 
be applicable to the defocusing case with a trapping potential, see also a related 
remark in |162j . 

It would also be interesting if one can find an asymptotic expansion, as e — ¥ 0, of 
the energy E\ of the first excited state (with least energy), as we did for the energy 
Eg of the ground state in Thcorcm ll.il The difference Ei—E g is of importance since 
it represents the "excitation energy" required to reach the first excited state from the 
ground state; it thus determines in some sense the stability of the ground state. (In 
the case of a convex bounded domain with Dirichlct boundary conditions, with the 
obvious modifications, this would provide evidence on the validity of a "nonlinear 
fundamental gap conjecture", see [24], [31]). 

If in addition the potential W is assumed to be even with respect to the coordinate 
axis, we have observed that one can construct a sign changing solution of (|1.22p _ . 
whose nodal set is the union of the coordinate axis, using the following strategy: 
Firstly, by minimizing the functional G-, described in Q1.23[) . over r\ £ W ' 
such that Wrj 2 € L^Rfj.), where R 2 + = {y = (y 1; y 2 ) : y x > 0, y 2 > 0}, for small 
e > 0, we obtain a positive solution in M.+ of the equation in (ll.22[) _ which is zero on 
the coordinate axis and approaches zero as |y| —¥ oo (we can see that the minimizer 
is nontrivial, if e is small, by adapting Example 5.11 in [3T] or Lemma 2.1 in [52"]). 
A solution u 2 defined in the entire space is then obtained using odd reflections 
through the lines y 1 = and y 2 = 0. The function u 2 is a solution of (jl.22|) _. 
whose 0-level set is the union of the two axis. Our construction parallels that of the 
well known saddle solution of the Allen-Cahn equation (jl.49j) . see [74]. The problem 
of existence and qualitative properties of saddle type solutions for the Allen-Cahn 
equation (not necessarily in two dimensions) has received a considerable amount of 
attention in recent years, see [49], [50], [51], [144] . We wonder if an analogous study 
can be conducted for the saddle type solutions of (|1.22p _ that we just described. 
Can one rigorously verify the formal prediction that 

u 2 -> sign{y 1 y 2 } v /(A - W)+ , 

say in L 2 (E 2 ), as e — > 0? The finer structure at the junction points on the axis, 
where W = A, may be demonstrated by a solution v of the following elliptic problem: 

v xx + Vzz - (x + v 2 )v = 0, x e M, z > 0, 
< v — \/ — x — > as x — > — oo; v — > as x — > oo, (6-3) 

v = if z = 0; v — V(x) — > as z — > oo, 

where V denotes the Hastings-McLeod solution as usual, which seems to be of 
independent interest. The above can be generalized to the case of arbitrary even 
space dimensions. Let us also note that, if the potential trap W is radial and two- 
dimensional, our construction can easily be generalized to obtain solutions u^. of 
(ll.22l) _ with N = 2, for small e > 0, whose zero level set has the symmetry of a 
regular 2fc-polygon and consists of k straight lines passing through the origin (see 
[15] for the corresponding solutions of (|1.49[) ). 

In the case where the potential W is, say, two-dimensional and symmetric with 
respect to the coordinate axis (as in the above paragraph) but the equation of 
(|1.22p _ is posed in R 3 , motivated from a definition in |102) . we can also consider 
"tick" saddle solutions: As before, minimizing the functional G- over rj G W ' (f2) 
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such that Wtf € L l (Cl), where fi = {y 1 > 0, y 2 > 0, -D < y 3 < £>)}, £> > 0, 
for £ < e(-D) sufhciently small, yields a positive solution of the equation which 
vanishes on dft = {y 1 = 0, y 2 = 0, y 3 = ±D}. By odd reflection with respect 
t° {Yi — 0j v 2 — 0) e (~^,D)}, and then with respect to the planes y 3 = 
(2k + l)D, k g Z, that solution can be extended to the whole of R 3 , yielding an 
entire solution of the equation in (| 1 . 22|) _ which has a saddle structure on each plane 
y 3 = constant and is periodic of period AD in the y 3 variable. 

Very recently, del Pino, Musso, and Pacard [53] studied entire solutions of the 
Allen-Cahn equation (|1.49l) which are defined in 3-dimensional Euclidean space and 
which are invariant under screw-motion. In particular, their nodal set is a helico'id 
of R 3 . We believe that, for sufficiently small e > 0, similar solutions exist for the 

3-dimensional defocusing Gross-Pitaevskii equation in (|1 .22|) with a 2-dimensional 

radial potential W (A, W satisfying our usual assumptions). What is the asymptotic 
behavior of these solutions as e — » 0? Can some results of [82] be generalized in our 
context? 

We believe that, if the potential W is restricted to the radial class, the approach 
of the current paper can also be applied to the study of the e — ¥ limiting behavior 

of vortex solutions of the NLS equation (| 1 . 2 1 [) . see [142] or [185] . namely solutions 

of the form 



where (r, 9) denote the polar coordinates in R 2 . Hopefully, the obtained estimates 
could be used to prove the, indicated by numerical evidence [137] (for the case of 
the model harmonic potential), orbital stability of u\ in the time evolution of the 
Gross-Pitaevskii equation, and thus answering the question raised in the end of the 
recent paper [175] . 

Non-degeneracy conditions of the form (|1.7j) are common in the study of transi- 
tion layered solutions of elliptic equations with bistable nonlinearity, see [95] , [201] . 
In that context, the surface <92?o represents the interface of the layer. It turns out 
that, in some cases, the aforementioned conditions can be removed completely (see 
[70) . [77]). In particular, the interface may be non-smooth or intersect the boundary 
of the domain. Motivated from this, we believe that one can show that rj e — > V 'A + 
uniformly in M 2 , or at least in compact subsets of M. 2 \dT>Q, as e — > 0, without as- 
suming condition (jl.7l) . (Here rj £ denotes the minimizer of G e or the ground state of 
(ll.21l) _). In this regard, we refer to [55] Prop. 3.16] for a related result (for (jl.22l) _ 
with q = 2). 

Is there a 'T-Convergence" theory [141] for f|l . 1[) . relating local minimizers of the 
limit functional (|1.14p to local minimizers of Ijl.ljl . as e — > 0? 

We wonder if, besides the one-dimensional profile V(x), there is a (genuine) two- 
dimensional one v(x, z) that could be used in (|3.10l) . In view of (I3.1ip . f|3.13[) . (|3.14p . 
and the matching conditions with y/ A + , the profile v should satisfy 

f P~ 2 v zz + v xx - v{v 2 + x) = 0, (x,z)eM. 2 , 




n = ±1,±2, 



(6.4) 




v — > as x — > oo, 



with v being £ £ /e s -periodic in z. As in [75] (see also [135] ). after a simple transfor- 
mation of the z independent variable (the x variable remains unchanged), abusing 
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notation, we are led to the problem: 

v zz + v xx - v(v 2 + x) = 0, (x,z)eR 2 , (6.5) 

v — y/—x — > as x — > — oo, v — > as x — > oo, (6-6) 

with v being £ e /si -periodic in z, where £ e = f " f3 e (9)d9. Uniqueness of positive 
solutions to the above problem does not seem to follow from the approach of Brczis 
and Oswald, as in [125] . since solutions are unbounded (compare with f|6.3[) ^1 . In fact, 
we believe that entire solutions of the equation in (|6.5j) should satisfy the growth 
estimate v(x,z) = 0(|x|2). Moreover, it is not clear how to adapt the uniqueness 
result of Brezis [42 . On the other hand, motivated from (|1.30l) . it is natural to seek 
solutions such that 

v x <0, (x,z)eR 2 . (6.7) 
It is irresistible to compare problem (|6.5|) . (|6.7I) with the famous De Giorgi conjec- 
ture for the Allen-Cahn equation (jl.49l) . see for instance [87], [107] . In this regard, 
it is interesting to investigate whether there are genuine two-dimensional solutions 
v of problem (|6.5|) . (|6.7[) or not. Note that solutions of the latter problem could be 
unbounded and, in particular, so could be v z (see |36[ 137] where a similar difficulty 
arises). We point out that the space dimension usually plays a very important role 
in these type of problems. Another direction could be to investigate the same ques- 
tion for stable solutions of (|6.5p . in the sense of (|3.75[) . see also a remark in pg. 79 
of the review article [72 . Actually, using the method of [87] , one can show that any 
solution of (|6.5|) . (|6.7[) is stable. A variation of these questions could be to consider 
problem (|6.5p - (|6.6[) , with the asymptotic behavior in (I6.6[) being uniform in z G M, 
along the lines of the so called Gibbons conjecture (see also [55]). We remark that 
in this case, as in |107j . the method of moving planes [108 can be applied to show 
that dnHJl holds. 

Many recent papers deal with the study of semiclassical ground states for the 
focusing (|1.22|) + . where the potential W(y) — A is positive but decays to zero, as 
|y| — > oo, at most like |y| -2 (see for instance [2 07) and the references therein). Can 
one study the defocusing case under analogous conditions? For e fixed, a related 
existence result may be found in [3j. 

Suppose, for simplicity purposes, that W is an even, double- well, one-dimensional 
potential (for example as in [124] or |139j ). say W(y) = (y 2 — l) 2 . What happens 
in the "degenerate case" when A equals the local maximum of Wl Assuming that 
W"(0) < 0, suitably blowing up at the origin, we expect that the fine behavior 
of solutions, as e — » 0, near the origin should be determined by a solution of the 
problem: 

r v " = v(v 2 + K^x 2 ) = 0, iel, 

Notice the similarities of the above problem with (|B.2[) below. Note also that in 
the case of a symmetric double-well potential (for any A) formulas (|6.2p do not hold 
due to tunneling effects, see for instance [104] . |121j . 

Finally, we believe that similar studies can be conducted in the case of the "ex- 
terior" problem, where A > W outside of a bounded domain and A < W in its 
interior. It is natural to assume that A — W — ^c>0as|y|— ^ oo, and consider 
the Gross-Pitaevskii equation (|1.2ip _ with boundary conditions \u(y, t)| — > c~ 
as |y| — > oo. The approach of [33] does not yield uniqueness of bounded ground 
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states, namely solutions of the equation in (| 1 . 22[) coupled with the aforementioned 
boundary conditions, if N > 3 (compare with Remark 13.1 2 j) . and one has to apply 
a sophisticated "squeezing" argument (see [84]). Let us mention that the stability 
of standing wave solutions of Gross-Pitaevskii equations, considered with nonzero 
boundary conditions at infinity, is a very active field of current research, see for 
instance [40] . 

Remark 6.1. It is worthwhile to mention that ifO < W(y) — A < c(l + |2/|) 2+d , V € 
M. , N > 3, for some positive constants c,d, then (jl.22p _ has infinitely many 
bounded solutions with positive lower bounds (see [136) ). 

Appendix A. A-priori estimates for the linearized operator based on 

THE NON-DEGENERACY OF THE INNER PROFILE 

Here we will provide an alternative, more natural, proof of the important Proposi- 
tion l3.5l that does not require knowledge of lower bound (|3.2ip , whose proof is rather 
technical (recall Remark |3.7|) . but instead relies merely on the non-degeneracy of the 
Hastings-McLeod solution V. This proof has the flexibility to deal with problems 
where the corresponding inner profile V is non-degenerate but the corresponding 
lower bound (|3.2ip may be hard to establish or fails (see Remark IB .41 below for an 
example where the latter case occurs). The latter situation certainly occurs when 
trying to construct unstable solutions (with respect to the parabolic dynamics) in 
related problems, see [13411135) . Let us also point out that it is not clear to us how 
to generalize the last part of the proof of (j3.21[) in [104) to the case of arbitrary 
power nonlinearity, as in Proposition IB . 1 1 below (see also Remark IB . 1 1 below) . 

PROOF OF PROPOSITION [331 Observe that it suffices to show the following a- 
priori estimate: There exists a constant C such that if e is sufficiently small, ip <E 
X n C 2+a (M. 2 ) , and / e X n C Q (M 2 ), < a < 1, satisfy 

« = /, (A.l) 

then 

IMU«(R») <C||/IU»(M=)- (A.2) 

To this end, as in the one-dimensional related problem treated in [191) Prop. 5.2], 
we will argue by contradiction. We remark that this indirect method has been 
used extensively in the study of elliptic singular perturbation problems involving 
transition and spike layers, see [2"U] . 

Firstly note that without knowledge of the validity of (|3.21[) . relation (j3.56[) 
would be 

!ces|x|, if L<|a:|<fe~§, 

(A.3) 
c + c|e3 y\ p , otherwise, 

for small e > 0, having increased the value of L if necessary. 

Suppose now that there exist sequences e n — > 0, ip n e X n C 2+Q (R 2 ), /„ e 
X n C Q (R 2 ) such that 

C(ip n ) = Aip n - e„ 3 (iu 2 ap - a(e%y)) <p n = /„ 

V ' (A.4) 

||<Pn||z«(R2) = 1 and ||/„||i=o (K 2) -> 0. 
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Keeping in mind that ip n — > as \y\ — > oo, we may assume that there exist y n G R 2 
such that, without loss of generality, we have 

tPniVn) = 1, Vyn(yn) = 0, A(p n (y n ) < 0, n > 1. (A.5) 
From (|A.3[i (|A.4|) . we obtain that 

dist(y„,f £ J <C, n>l, (A.6) 
for some (generic) constant C independent of n > 1. Thus, abusing notation, we 
/3 _1 (£^z„)a;„, z„j with |x„| < C, z„ € [0,£„ 3 4„)- Therefore, 
passing to a subsequence, we may assume that 

2 

a: n — > and e£z„ — » e [0, 4)]- (A. 7) 

Recalling (13. 6[) . (13.101) . in terms of coordinates (x, z), the equation in (IA.4|) takes 
the form 

2 „ _ 2 / 2 \ 

(<£n)** + ^(^X^rOxx + Blfon) ~ £n 3 ( 3u 2 p - y) J fn = fn, (A. 8) 

in the neighborhood of the curve f Sn described by j |x| < (5 e„ 3 , 2: € 0, e n 3 ^ | , 
where Si is the differential operator: 

Bifa) = eip"p- l xtp m + el(/3') 2 /3- 2 a;V^ + 2e^P , j3- 1 xip mj , + B^ip), (A.9) 



and B\ is the differential operator in (|3.7I) where derivatives are expressed in terms 
of formulas (|3.11[) and s replaced by fi~ x x. By (|3.18|) . the first relation in (|3.40|) . 
and working as in ((3T57J) , we obtain that 

3u 2 ap - a{ely) = elp 2 {elz) (3V 2 (x) + x) + 0{el){x 2 + 1) (A.10) 

uniformly in the region described below (|A.8|) . as n — > 00. Making use of (|A.4[) - 
(|A.10p . and a standard compactness argument, passing to a subsequence, we find 
that 

y n ^<f in Cf oc (R 2 ), 

where <j> satisfies 

0«+ $ ffi(z*)(W 2 (x)+x)0 = O 1 (x,z)eR 2 , (A.ll) 
||0IU«(B») = 1 (<Kx m ,z.) = l). (A.12) 

Since 

3V 2 (x) + x -> 00 as a; -> ±00, (A. 13) 

a standard barrier argument, as in Q3.72p . and elliptic estimates (109) . yield that 
there exists a constant C such that 

\V(f>(x, z)\ + \<f>(x, z)\ < Ce~W, (x, z) G M 2 , (A.14) 

(see also Lemma 7.3 in [801 ) . Let (/!!, -0!) denote the principal eigenvalue-eigenfunction 
pair of 

-M(i^) = -il)"+(3V 2 (x)+x)^ = fj,ip, V(±oo) = 0. (A.15) 

Without loss of generality, we may assume that ^\ is positive. Furthermore, we 
have 

Hi > 0, (A.16) 
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as testing against V x < readily shows (see |191j ). (By virtue of (|A.13|) and 
Theorem 10.7 in [121] . the spectrum of — Ai, in L 2 (R), consists of simple eigenvalues 
/Ui < A*2 < • • • with fj,i — > oo). Now let 

/oo 
<f>(x, z)ipi (x)dx, 
-oo 

where <j) is as in (|A.11|) . From (|A.11|) , (|A.14|) . and (|A.15I) with ij) = ip\, fi = (J,i, wc 
calculate that 

= 0o(**) Ho [-^ + ( 3 ^ 2 (*) + *) 0] <M^x 
= /3 2 (z*) [-^' + (W 2 (x) + x) ^] 

and 

|$| < C, zel. 

From (|A.16|) . and the above two relations, it follows at once that $ is identically zero, 
which contradicts the previous relation ([A. 121) . Consequently, we have established 
the validity of the desired a-priori estimate (|A.2I) . 
The proof of the proposition is complete. □ 

Remark A.l. By adapting Lemma 5.3 of [191] . we can show that relation iA.l]) 
implies that || V 3 !! (M 2 ) < C||/IU 2 (R 2 ) f or some constant that is independent of tp, f 
and small e > 0. In fact, as in [101] Thm. 1.2], we expect that more general a-priori 
estimates of the form \\(p\\lp(r 2 ) < Cp,q £Qp ' <! ll/l|L < !(R 2 ) hold true. 

Appendix B. Around the Hastings-McLeod solution of the 

PAINLEVE-II EQUATION 

In this appendix we will provide a new proof of the existence of the Hastings- 
McLeod solution V of the Painleve-II equation p. 151) . Moreover, we will establish 
various qualitative properties of the solution that are required for the singular per- 
turbation analysis. In contrast to the original proof of Hastings and Mcleod [113] 
(see also [116] ). where a shooting argument was employed, here we will use an upper 
and lower solution argument, which in principle is not restricted to ODE problems. 
Even though such an approach was successfully applied to this problem recently 
in |18j . and very recently in [202} , in our opinion our construction is more flexible 
and intuitive. The main advantage of our proof, compared to those of the latter 
references, is that, in the process, we also establish existence and uniqueness of a so- 
lution of problem (|3.107]h which seems to be a new and useful result (recall Remark 
13.13]) . Although a sizable literature has been devoted to the study of the Painleve 
equation (see [9U1 11291 1149] ). we understand that the solution of this problem was 
not previously known. 

Our choice of lower-solution is motivated from [71] where, in particular, the 
authors treat the problem 

u" = \u\ p — x, x > 0, u(0) = 0, x~pu(x) — y 1 as x — y oo , (B.l) 

where p > 1 (recall the discussion in the third part of Subsection II .2]) . see also [116] 
Sec. 3.2] and [122] for the case p = 2 which is the Painleve-I equation. On the other 
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hand, our choice of upper-solution is motivated from [190] where, in particular, the 
authors treat the problem 

u" = u 2 - x 2 , x € M, u(x) - \x\ -> as \x\ -> oo, (B.2) 

see also Remark IB .81 below. 

The notation in this appendix is independent of the rest of the paper. 

As in |113j . see also [32] pg. 200], we will prove the following more general result: 

Proposition B.l. Given p > 1, there exists a unique nonnegative solution U of 

- u" -xu+ \u\ p u = 0, x E R, (B.3) 

such that 

u(x) — > as x — > — oo; u(x) — x~p — > as x — > oo. (B.4) 
Furthermore, we have that U' > in R, and 

U(x) = O (\x\~^e~i^\ as x -> -co; U(x) - x* = O ^a;^ 2 ^ as x ->• oo. 

(B.5) 

The solution U is non-degenerate in the sense that there are no nontrivial bounded 
solutions of the problem 

<t>" - \{p + l)U p - x] 4> = 0, lei 

Note that (|1 .37|) falls in the above class of problems by means of the transforma- 
tion x — ¥ —x. 

Remark B.l. The results of this appendix can be used in extending the results of 
the current paper, and treat the defocusing (|1.22j) _ with arbitrary nonlinearity ex- 
ponent q > 2. In particular, the considered model of the latter equation, with non- 
linearity exponents 7/3 < q < 3 (in one space dimension and the model harmonic 
potential), is particularly relevant to the physics of BEC-BCS (Bardeen-Cooper- 
Schrieffer) transition in ultracold Fermi gases (see [206] ). 

As a stepping stone towards the proof of the above proposition, we will first prove 
the following result which, as we have already mentioned, is of interest in its own 
right. In particular, the solution U + below will form the basis for our construction 
of a lower-solution to problem (|B.3l) - (|B.4p . 

Proposition B.2. Given p > I, there exists a unique solution U + of the problem 

' -u" = xu- U p +1 , x>0, 

< «(0) = 0, (B.6) 

< u(x) < xp , x > 0. 
Furthermore, we have that U^(x) > 0, x > 0, and 

1 / l-3p \ 

U+(x) — xp — O \x p j as x — >■ oo. (B.7) 

Proof. It is easy to check that xp is an upper-solution of (|B.6[) . while 5x[k,k+it] sin(x— 
K) is a (weak) lower-solution provided that K > 2 and < <5 < 1 (here x denotes 
the characteristic function). (We refer the reader to [53] for more information on 
piecewise smooth weak upper/lower-solutions). From now on we fix such a K, say 
K = 7r. Then, by a well known theorem [53], for every < S < 1, there exist 
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solutions m 2 of (|B.6|) such that Sxw^-k] sin(x — 7r) < u\ < m 2 < , with the 
property that any solution of (IB.6|) such that 8x{k,2tt\ sin(x — 7r) < u < xp , satishes 
Mi < m < m 2 . (Note that U\, m 2 depend on S). 

For any nontrivial solution of (|B.6|) we have u" = u(u p — x) < 0, x > 0. Thus 
it' is non-increasing. As a result, u' 4 a as x -> oo. Here a may be minus infinity 
We claim that a — 0. In fact, if a < 0, then m(:e) < for x large, which is a 
contradiction to u > 0. If a > 0, then u(x) > |i for i > large, which is a 

contradiction to m(x) < x~p. Thus a = 0. Consequently m' > 0. Actually, since u' 
cannot be constant over a nontrivial interval, we find that u' is decreasing and thus 
v! > 0. 

In order to proceed further, we will show that problem (IB. 61) has a unique so- 
lution. Inspired by the work by H. Brezis and L. Oswald in [33] (see also [1711 
Sec. 8.5.2]), we will study the quotient of two solutions (for a different approach 
see Remark IB . 21 below) . Suppose that (|B.6[) has two solutions u\, m 2 . Then, we can 
find < 8 < 1 small such that 8x[n,2-n\ sin(x — 7f) < ttj < X?, i = 1, 2. Thus, from 
the previous discussion, we infer that u\ <Ui < 112, % = 1, 2. We only need to prove 
that Mi = U2- From (|B.6|) . we find that 

<> < 
u 2 ~ Ul 

The above inequality implies that the function u' 2 Ui — U2u[ is non-decreasing, and 
so 

(u' 2 Ui — U2U l 1 )(x) > (u 2 Ui — U2lt'i)(0) =0, x > 0, 

which in turn implies that the function — is non-decreasing. Therefore, we get 



^M. > Km ^aM = x > 0, 

Mi (a;) K-i-o+ Mi(x) u'i(0)' 

by L'Hospital's rule (recall that m'^O) > 0). From M2 > Mi and mi(0) = m 2 (0) = 0, 
we know that m' 2 (0) > m^O). Suppose that m' 2 (0) > m^O). We have 

ui(0) mUO) i 

M 2 (0) M 2 (0) 

So 

U 1 = Ul(u1 — X) < Ml ' ' 



f 



X < 



for x > large, and some constant c > 0, since m'^O) < m 2 (0) and m' x > 0. It follows 
that u'± < for x > large. This is a contradiction to m'^it) > 0, x > 0. We 
conclude that 1^(0) = m' 2 (0), which gives mi — m 2 . 

Let m denote the unique solution of (|B.6|) . Adapting an argument from [71j . we 
will show that 

x~pu(x) — > 1 as x — > 00. (B.8) 

From m'(oo) = 0, we obtain 

/>oo 

u(x) (x - vP{x)) dx = - u"(x)dx = u'(0). 
Jo 

Hence, we can choose Xi — > 00 such that m(x;) (a;^ — u p (xi)) — > as i — >• 00, and 
recalling that m' > 0, we find that 

Xi — u p (xi) — > as i — > 00. (B.9) 
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Now, for any small 9 > 0, we claim that 

u(x) > (1 — 9)x~p , 

for x > large. Suppose that the claim is not true. Then, there are < 9 < 1 and 
Sti — > oo such that 

u(xi) < (1 - 9)xf . (B.10) 

It is easy to check that there is a T > large such that, for x > T, we have 

/ i \ " i / i \ P+ 1 

-((1-9)xp) < x(l - 9)xp - ( (1 - 6)xp J . (B.ll) 

By (|B.9|) . we can choose a constant f > T, such that w(T) > (1 - 0)Tp. Define 
v (x) = u(x) if x £ [0, T]; v(x) = max fu(o;), (1 — 9)x~p^j if x G [T, oo). Then w is 

continuous, < v < x?, and is a (weak) lower-solution of (|B.6[) . in view of (IB.11|) 
and [33]. As a result, from [1681 Thm. 2.10], problem (|B.6[) has a solution u* with 
v < u* < xp. On the other hand, since v > u, and u ^ v (by (|B.10[) ). we find that 
u* 7^ u. This contradicts the uniqueness of the solutions of (|B.6[) . 
It remains to show the validity of estimate ()B.7|) . We have 

(uP — x ) l l 

u" = u± t-{u -xp) = q{x)(u -xp), x>0, (B.12) 



with 

q(x) 



pasi-> oo, (B.13) 



by (|B.8p . (Note that, by the maximum principle, we get u{x) < X? , x > 0). Let 

i 

w = u — xp , x > 0. 
Then, from (|B.12[) . we obtain that 



"M f \ 1 f 1 

w + q(x)w — 

P \P 



- 1 x~p =0, x > 0. (B.14) 



We claim that there exist L, M > sufficiently large such that, for every S £ (0, 1), 
the function 

— —3 — 

w_(x) = —Mxp — 5xp , 

is a lower-solution of (|B.14[) in [L, oo). Indeed, thanks to (|B.13|) . for every S £ (0, 1), 
we find that 



+ q[x)w- i (± - 1J xl' 



M (I - 3] (| - 4) xp- 5 + 51 (I - l) xi- 2 - fAfxi- 2 - |«5x? +1 - I (l - 1 ) , ~ 



Mi?" 5 



<0, 



i-3jU-4j + (i-f). 

and L is Si 

4 
p 

w{L) < w(L) for every S £ (0, 1). In view of (|B.8|) . for every 8 £ (0, 1), we have 
w(x) — w{x) — —Mx~p~ 3 — Sxp — (u — xp) < xp [—5 — (x~pu — — > — oo 



for x > L, provided M > - and L is sufficiently large (independently of M). We 
fix such an L > 0, and choose M > - large, depending only on L, such that 
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as x — > co. Thus, by the maximum principle, we deduce that w(x) < w(x), x > L, 
i.e., 

— a;p > —Mxp ' — Sx" , x > L, < <5 < 1. 

By letting S — > 0, we deduce that tt(ai) - if > —Mxv , x > L. The validity of 
(|B.7|) follows at once from this lower bound and the upper bound u < x? . 

The proof of the proposition is complete. □ 

Remark B.2. An alternative way to establish uniqueness for W.6\) is the following: 
Suppose that V+ solves IB. 0(1 . Then, it is easy to see that XV+, A > 1, is a family 
of upper solutions of KB.6\) such that AV+(0) = U+(0) = 0, and W+ — U+ — > co as 
x — y co if A > 1. Moreover, since V, (0) > 0, we have that W+ > U+ in (0, co) 
iJA > 1. Therefore, by Serrin's sweeping technique (see |184[ pg. 40],), we get 
V+ > U+ in [0,oo). Similarly, we can show that V+ < U+ in [0, co). Consequently, 
we get that V+ = U+ . 

Remark B.3. In relation with the problems mentioned in the third part of Sub- 
section \ 1.021 it would be of interest to generalize Proposition \B.2\ in the following 
direction: Study solutions u : M 2 — > C (if they exist) of the problem 

Au + (\y\ - \u\p)u = 0, yeR 2 , 

u(0)=0; \y\~ »\u\ I as \y\ oo. 

We cannot resist to compare the above problem with the well known 
Au+ (1 - \u\ 2 ) u = 0, yeM 2 , 

u(0) = 0; |u| -+1 as \y\ -> co, 

see [1 20] . |170j . and the references in the research monographs [39], [171] . [183 . 

The following proposition plays an important role in relation with Remark 13. 131 

Proposition B.3. The solution U+ of \B. 6\) is non-degenerate in the sense that 
there are no nontrivial bounded solutions of the problem 

<p" - [0+ l)U p + -x] <j) = 0, x > 0, </>(0) = 0. (B.15) 

Proof. Suppose that there exists a nontrivial bounded solution (j> of (|B.15p . The fact 
that (p+l)?7i — x — y co as x — y co , easily implies that <fi has a finite number of zeros 
in [0,co), all of them simple, and that |</>|, \(j>'\ decay to zero super-exponentially 
as x — y co. Let ro > be the largest zero of cf>. Without loss of generality, we may 
assume that 4>'(r ) > 0. Differentiating (|B.6I) . multiplying the resulting identity 
by </>, then multiplying ([B.15[) by U' + , subtracting, and integrating by parts over 
(ro,co), we readily arrive at 

/•oo 

U' + (r )<P'(r ) = - / U + 4,dx. 

Jr 

However, this is a contradiction to U' + > 0, <f>'(ro) > 0, and 4>{x) > 0, x > ro. 
The proof of the proposition is complete. □ 

Remark B.4. Note that, since U+(0) — and p > 1, the potential (p + 1)J7? — x 
of the linear operator in the righthand side of \B.15\) is negative for small x > 0. 



THE GROUND STATE OF A GROSS-PITAEVSKII ENERGY 



(51 



We can now give the 

PROOF OF PROPOSITION IbTTI Let U+ be as in Proposition IB. 21 it is clear that 



u(x) = { (B.16) 




is a (weak) lower-solution of (IB.3|) - (IB.4I) . 

Next we will construct an upper-solution of (IB.3[) - (|B.4[) . for a different construc- 
tion we refer the interested reader to Remark IB. 101 below. Let 




u (x) = { (B.17) 



Then, fix a continuous function <fi > such that uo + 4> € C 2 (R) and <p(x) = if 
\x\ > 1. Let ji\ > 0, tpi > denote the principal eigenvalue and the corresponding 
i^-normalized eigenfunction of 

-f + ((p + lK-x)^ = ^, ^ei 2 (R). 

Such fii, ipi exist, since the potential 

( px, x > 0, 

Q(x) = (p+l)u p -x = I (B.18) 

[ —x, x < 0, 

clearly satisfies inf xG R Q(x) < ]im.ini x ->.± 00 Q(x), see |178) (see also |121l Thm. 
10.7]). Furthermore, ipi, \ip"\ decay to zero super-exponentially as |x| — > +oo. 

More precisely, there exist constants c± and X- < 0, x + > such that 

ex P { ~ Jx± y/Q(t) - ^ dt } 
ipi(x) ~ c± 1 — as a; — > ±oo, (B.19) 



[QW-mi] 3 

see [521 Chap. IV, Thm. 14] and [THl Appx. A]. We can now define our upper 
solution for (|B.3|) - (IB.4|) as 

u = u + (f) + MQ D ipi 

with M, D > 1 large constants to be chosen, and Cd(x) — z(x — D), where z € 
C°°(M) is such that z = 1, t < 0; 2 = 0, t > 1; z' < 0, te (0,1) . We proceed in 
verifying that u is indeed an upper solution. In [0, 1], u = (uq + 4>) + Mtp\ and 



xu + w 



p+1 > -CM + cM p+1 > 0, 



provided M > is sufficiently large (C, c > are independent of large M). We fix 
such an M > 0. In [1, D], we have u = uq + Mipi and 

- u" - xu + u p+1 = -u'^-M^-xu a -Mxi' 1 + (u Q + M^ 1 ) p+1 

> -Mtp" - xu - Mxrpi + u p +1 + (p+ l)u^Mip! = > 0. 
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In [D, D + 1], we have u = u + MQu^i and 

- u" - xu + u p+1 = -u'^- M(Cd^i)" -xuq- Mx( D tp 1 + (u + MC D ij 1 ) p+1 
> -- ( - - l) x^- 2 - Ce~ x - xu + < +1 



P \P 

= -- (- - l)xp- 2 -Cer x > 0, (B.20) 
P \P J 

provided D is chosen large (C > is independent of D). In [D + l,+oo), we 
plainly have u — xp. Analogous calculations also hold in (— oo,0]. Consequently, 
the function u is an upper-solution of (IB.3|) - (1B.4I) . 

It follows from [1681 Thm. 2.10] as before that there exists a solution of (|B.3I) 
such that u < u < u. The second estimate in (|B.5j) follows at once from (|B.7|) and 
(|B.16|) : the first one follows from the fact that, for every solution of (|B.3|) that tends 
to zero as x — > — oo, there exists some constant c > such that 

u(x) ~ cAi(— x) as x — > — oo, (B.21) 

(recall the discussion leading to (JX2SI), and (|3.24p . 

We will show that v! > 0. We follow |18j . Since < u < x? for x > D + 1, as in 
the proof of Proposition [BT2] we obtain that v! — > as x — > oo. Moreover, it is easy 
to show that u' — > super-exponentially as x — oo. Since it is strictly positive 
(by the maximum principle), in view of (|B.3I) . we can write 



i.e., 

v vv — pu p v = —u, where v = u' . (B.22) 

u 

Since v — > as |a;| — > oo, it follows that if v is not strictly positive, then there exists 
Xq such that 

v(x )<0, v'(x )=0, v"(x o )>0, 
but this is impossible because (|B.22[) implies 

v"(x ) ~ pu p {x )v{x ) = -u(x ) < 0, 

and therefore we conclude that u' is strictly positive. The same conclusion can 
also be derived by adapting an argument from [202j . i.e., applying the maximum 
principle to the function — (keep in mind the second identity in (|3.29jl ). A more 
PDE approach is to apply the moving plane method |108j . starting from — oo, as in 

Uniqueness (of nonnegative solutions) for the problem (|B.3[) - (|B.4[) can be estab- 
lished in a similar manner as we did in Proposition lB.2l for the problem (IB.6|) : Again 
we suppose that there exist two distinct non- negative solutions m, u 2 of (|B.3|) — (|B.4|) . 
By the strong maximum principle, we deduce that both are strictly positive. Hence, 
there is some small 6 > such that Ui{x) > Sx[-K,2Tr] sin(x — 7r), i = 1, 2, for every 
x £ M. Recall that the function in the righthand side is a lower-solution of (|B.3|) ~ 
(|B.4|) . Moreover, both Ui, i = 1,2, will eventually lie below the graph of xp (note 
that u — xp is strictly convex as long as it is nonnegative). Thus, by virtue of (|3.24p , 
(|B.18|) . (|B. 191) and (|B.21|) . we can choose sufficiently large numbers D, M such that 
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Ui(x) < u(x), i = 1, 2, for every i?R. (This is a fine point that was not present in 
the uniqueness proof for (|B.6p ). Recall that the function in the righthand side is an 
upper-solution of (|B.3|) - (|B.4j) . Consequently, we may assume that u\{x) < U2(x), 
x € K, and it is easy to see that ^ is non-decreasing in R. On the other hand, it 
follows from the second relation in (|B.4[) that ^ — > 1 as x — > oo. So, we get that 
U2 < Mi in R which is a contradiction. 

Finally, the non-degeneracy of U can be derived as in Proposition IB.3I for the 
non-degeneracy of Z7+. 

The proof of the proposition is complete. □ 

Remark B.5. The fact that problem liB.3\ )~ fB.4\ ) has a unique solution, which is 
a stronger result, has been proven in |113) . 

In contrast, problem iB.l]) withp — 2, has exactly two solutions. Existence of two 
solutions has been established by Holmes and Spence |122) by a shooting argument 
(and in [71] for any p > 1, via the method of upper/ lower solutions and variational 
arguments, perhaps unaware of [122 ), where the authors also conjectured that these 
solutions were indeed the only ones. Their conjecture was settled, to the affirmative, 
by Hastings and Troy [114] . However, their proof was, as we discover now ( almost 
25 years later! ), much more complicated than necessary, and relied on some four 
decimal point numerical calculations. Motivated by an idea of ours from j!35j . where 
problem \B. Sty was shown to have exactly two solutions, we can give a truly simple 
proof of the uniqueness result of |114] as follows. We know from |71j . [122j that 
problem \B. 1\) with p = 2 has a unique increasing solution 14+ . Let U be any other 
solution, and let rj — IA+ — U . By an easy calculation, and the maximum principle, 
we find that rj has to be a positive solution of 

n" - 2U + {x)ri + rf = 0, x > 0, r/(0) = 0, rj(x) -> as x -> oo. (B.23) 

The key observation now is that the solution rj furnishes an odd standing wave 
solution of a focusing NLS equation of the form (jl.21|l +. with N = 1, q = 2. 
Thankfully, in the last years a lot of research and efforts of many authors have been 
devoted to the uniqueness of ground states of radially symmetric focusing nonlinear 
Schrodinger equations with non- decreasing potential (in r > 0) and power nonlin- 
earity, considered in the whole space, in a ball, or an annulus (see |48j for the state 
of the art). The problem of uniqueness of rj resembles more the case of the annu- 
lus and, having all those tools at our disposal which were not available at the time 
that |114) was written, we can infer that uniqueness as well as non-degeneracy of 
a positive solution rj of \B. 23\) follow directly from |911 Thm. 1.2]. (In the latter 
reference, it was assumed that the potential is strictly positive but it is easy to check 
that their proof works equally well for the case at hand, see also |48j . |132 ). The 
non- degeneracy property of the solution U, which follows readily, is a new result 
and, in the context of the original singular perturbation problem [197] (see also |71j . 
|73j ) is more useful than uniqueness (recall Remark \3.7\ ). 

To the best of our knowledge, the similarities between the singularly perturbed 
(multi- dimensional) elliptic problem in [71] , arising from the study of the Lazer- 
Mckenna conjecture, and the one- dimensional one in [114] . [122j . |197j . arising from 
the problem of vertical flow of an internally heated Boussinesq fluid with viscous 
dissipation and pressure work, were previously unknown. 

The proof of the Lazer-Mckenna conjecture in [71j consists of constructing so- 
lutions of the problem with arbitrary many (clustering) sharp downward spikes on 
top of a positive minimizer of the corresponding energy functional, as the small 
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parameter e > tends to zero. The aforementioned minimizer has a corner layer, 
along the boundary of the domain, whose profile is described by the positive solu- 
tion of lB~J\) . It is our hope that the techniques of the present paper, together with 
those already developed in [78j . |160] . can be used to construct new solutions, having 
corner layer profile described by the unstable solution of /IB. 1\) (is this unique and 
non-degenrate for p ^ 2?), at least when e stays away from certain critical values 
£\ ~> e 2 > ' ' ' > £j — ► (see [135] for a related problem). These solutions would 
be slightly negative in a small (e-dependent) neighborhood of the boundary, with 
their Morse index diverging as e — > (away from the critical numbers). Then, one 
could use the techniques of [85] . [200] in order to add downward spikes on top of that 
unstable solution. One may even be able to prove the existence of arbitrary many so- 
lutions, as e — > 0, which was the original assertion of the Lazer-Mckenna conjecture 
settled in [71] . just by the fact that the solution's Morse index diverges (recall also 
the discussion following i6.2\) ). This would constitute a proof of the Lazer-Mckenna 
conjecture, as treated in [71], [73], that is valid even for supercritical exponents. 

Remark B.6. The related boundary value problem 

u xx = u(u 2 — x), x > 0; u x (0) — 0, u — \fx — >• as x — > oo, (B.24) 

arises in the study of the superheating field attached to a semi-infinite supercon- 
ductor, for the construction of a family of approximate solutions of the Ginzburg- 
Landau system via the procedure of (formally) matching inner and outer solutions 
(see [571 1117] ). Existence for the problem \B.2J$ has been established by shooting 
arguments in [117] . and by topological ones in [172) . We can give a new proof of 
their results, valid for any power nonlinearity ( as in Proposition \B. 2\) . by slightly 
modifying the above proof of Proposition \B.1\ as follows: One still uses U+ as a 
lower solution of \B. 24\ ) (it is a solution and f/+(0) > 0, see [21],),' however in the 
construction of the upper solution we have to be careful to chose the function 4> 
such that (uq + cf>) x (0) < (and afterwards the principal eigenfunction subject to 
Neumann boundary conditions at x — 0J. It was shown in [112] that, without any 
assumptions at infinity, problem \B. 24\ ) has exactly one global positive solution. As 
a matter of fact, we expect that an analogous property holds true for the problem 
l(B.6\) , see also [42] for a related result concerning fB.2\) . 

Remark B.7. Since p > 1, the solutions that we have constructed in Proposi- 
tions \B.1[ \B.2\ have infinite energy ( more precisely, their derivative does not belong 
in L 2 (0,oo)). Nevertheless, we believe that one can also establish existence (and 
further characterize the solutions) for problems (B.S\) - fB.4\ ), /IB. 6\) . and \B. Wm by 
minimizing a suitable re-normalized energy functional, as in [170j . or minimizing 
the standard energy functional of iB.3\) in a large interval [-R, R] with boundary 
conditions u(—R) — 0; u(R) = R7 , and then letting R — > oo (see [37], [49], [82] for 
some related situations) . Similarly for the other cases. 

Remark B.8. If0<p< 1, then there exists a unique solution of \B. 3) - \B.l$ such 
that u > uo, where uo is as in HB.17\ ). We cannot resist to give a short proof of 
this, based on [134] . If0<p<l,itis easy to see that uq is a lower solution, while 
uq + ip, with if > defined by 



-tfixx + [(p + l)"o ~x](f= (u ) xx > 0, <^(±oo) = 0, 
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is an upper solution of KB.3\) - fB.4\ ). If p = 1, then uq is a weak lower solution, 
while uq + <j), with (f> > the unique continuous solution of 

-4>xx + [(p + - x] <f> = 0, 0,(0-) - 0,(0+) = 1, 0(±oo) = 0, 

is an upper solution of HB.3\) - fB.4\ ) (recall HB.18\) ). Existence of the desired solution 
follows at once. Uniqueness follows simply by taking the difference of the equations 
satisfied by two pairs of solutions. 

The case p = 1, which is equivalent to HB.2\) . has received considerable atten- 
tion lately, mainly since it appears in the study of the spatial segregation limit of 
competitive systems [62] , [127] (see also |123j . [190] ). Interestingly enough, we have 
noticed that problem /IB. 2]) also describes the corner layer profile of solutions in the 
paper [115] by Hastings and McLeod, in the case where the constant c therein, which 
arises from an integration, is chosen to be zero rather than strictly positive. 

Remark B.9. Similar results should also hold true for the equation 

u" + x\x\ s u- \u\ p u = 0, xeR, 

where p, s > 0. 

Remark B.10. As we have already seen in Remark \B.8\. there exists a unique 
solution U > max{0, x} of the problem 

u" = pu{u - x) = 0, (B.25) 

such that U — > as x — > — oo; U — x — > as x — > oo. (We have found out in |135j . 
by arguing as in Remark \B.5\ above, that there exists exactly one more solution U- 
of W.25\) which satisfies the boundary conditions, and in fact U- < max{0, x}). 
Actually, the solution U is the unique (global) solution of W.25\) such that u > 
^, x G K. This follows at once from the fact that the equation in KB.2\) has a 
unique nonnegative solution, see [42] . Moreover, since p > 1, it is easy to see that 
the function U~p is an upper- solution of HB.3\) - fB.4\ ). 
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